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ABSTRACT 


An overview of space electric propulsion (SEP) is presented. Methods of 
throttling the power levels of electrostatic and electromagnetic thrusters are 
discussed. Particular attention is given to the concept of thermally-throttling 
propellant flow using the temperature-viscosity characteristics of xenon gas. 
The thermoproperties of xenon gas as a function of temperature are 
determined, and the flow regimes of the propellant at the mass flow rates of 
interest are studied. The propellant flow is presented separately as Fanno flow 
and as Rayleigh flow, and then those combined effects are considered. A 
method for predicting the performance of thermally-throttled systems is 
presented. Uncertainties in modeling real-world thermal throttling systems 
are discussed. The possible use of thermal throttling characteristics as a means 


of propellant pressure regulation is also examined. 
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I. INTRODUCTION 


A. SPACE ELECTRIC PROPULSION 


Electric propulsion devices for space applications use electrical energy 
(from solar cells, batteries, fuel cells, etc.) as a means to heat or to eject 
propellant mass, resulting in thrust. While there are several methods of 
electric propulsion either in use or in development, all can be described as 


one of three fundamental types (Sutton, 1992): 


1. Electrothermal thrusters. Electrical energy is used to heat the 
propellant, which is then thermally expanded and accelerated to 
supersonic speeds through a nozzle, as in a typical chemical rocket. 

2. Electrostatic thrusters. Thrust is produced by accelerating electrically 
charged propellant particles through an electrostatic field. 

3. Electromagnetic thrusters. Electrical energy is used to produce a 
plasma (a high-temperature electrically neutral gas containing both 
electrons and positive ions) which interacts with electric and 


magnetic fields to accelerate the propellant. 


The three basic subsystems of a typical space electric propulsion system 
are (1) a power processing unit (PPU), which contains electrical interfaces 
between the propulsion system and the spacecraft, and sends power and 
command signals to the thruster and propellant flow control signals to the 


propellant subsystem; (2) a propellant management assembly (PMA), which 


provides a means of storing, throttling and delivering the propellant; and (3) 
one or more thrusters, which convert electrical energy to thrust through one 


of the means discussed above. 


T Electrothermal Thrusters 


There are two basic types of electrothermal thrusters currently in use. 
The first and simplest is known as a resistojet (Figure 1), which uses electrical 
current to directly heat metal components such as coiled wire, tubing or fins. 
These components then heat the propellant through radiation and 
convection as it flows over them. In the example shown below, an increase of 
about 28% in specific impulse over a standard catalytic thruster is achieved by 
electrically augmenting the enthalpy of the decomposition products of 
hydrazine using a high-power heater in a second-stage vortex chamber; thrust 
results when the higher energy propellant is accelerated through the exit 
nozzle. Resistojet performance is limited primarily by the capacity of the 
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Figure 1. Resistojet thruster cross-section (From Agrawal, 1986) 
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thruster's structural components to withstand high temperatures. This 
limitation is overcome by the second basic type of electrothermal thruster, 
known as an arcjet. These use an attached arc within the nozzle to directly 
heat the propellant stream to temperatures much higher than that of the 
thruster body. Arcjet thrusters are limited primarily by instabilities in the arc 
and erosion in the nozzle. 

Resistojets can produce thrust levels of 300-500 mN at a specific 
impulse of about 300 seconds and an efficiency of 65-9076. Arcjets can produce 
thrust levels of 200-700 mN at a specific impulse of 400-1600 seconds and an 
efficiency of 40-6575. Electrothermal thrusters have been flight-tested on a 
number of spacecraft, including Intelsat V. NASA currently has two 


electrothermal thrusters under development: a 600 s Isp, 2 kW arcjet for 


commercial communications satellite north-south stationkeeping, and a 500 s 


Isp low power arcjet thruster (LPAT). The U.S. Air Force currently has both 30 


KW and 1 KW arcjet thrusters under development. 


2. Electrostatic Thrusters 


Electrostatic thrusters use electrostatic fields to accelerate charged 
particles in a linear direction, producing thrust. While electrons are relatively 
easy to produce, their small mass makes them impractical for use as a 
propellant. Accordingly, electrostatic thrusters use heavy atoms charged as 
positive ions (since neutralization with electrons is practical) or charged 
colloids (liquid droplets). These thrusters can be categorized by their source of 


charged particles as follows (Sutton, 1992): 


1. Electron bombardment thrusters. Positive ions are produced by 
bombarding vaporized or gaseous propellant (such as xenon or 
mercury) with electrons emitted from a heated cathode (Figure 2). 

2. lon contact thrusters. Positive ions are produced by passing 
propellant vapor (usually cesium) through a hot contact ionizer 
(usually tungsten). 

3. Field emission (colloid) thrusters. Tiny liquid droplets of propellant 


are charged by passing them through a corona discharge. 
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Figure 2. Electron bombardment ion thruster (From Agrawal, 1986) 
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Electron bombardment ion thrusters can produce thrust levels of up to 
100 mN at a specific impulse of as much as 5000 seconds and an efficiency of 
50-65%. Ion contact thrusters can produce thrust levels of up to 500 mN at a 
specific impulse of as much as 5000 seconds and an efficiency as high as 90%. 
Field emission thrusters can produce thrust levels of up to 200 mWN at a 
specific impulse of 4000-6000 seconds and an efficiency of more than 95%. 
There are currently no U.S. development efforts directed toward either ion 
contact thrusters or field emission thrusters. Electron bombardment thrusters 
have, however, been flight-tested on a number of U.S. spacecraft, and there 
are several programs pursuing further development of these thrusters in the 
United States. Chief among these, NASA is currently developing a 30- 
centimeter xenon ion engine for the NSTAR (NASA Solar Electric 
Propulsion Technology Application Readiness) program, which will serve as 
the powerplant for the New Millennium spacecraft scheduled for launch 
beginning in 1998. The U.S. Air Force Office of Scientific Research (AFOSR) is 
also involved in electrostatic thruster development, funding research of a 
variant which would use C60 (Buckminsterfullerene) as a propellant. Heavy 
particles such as this are desirable in electrostatic thrusters, where the thrust 
per unit area increases as the square of the particle mass to charge ratio; such 
particles contribute the desirable operating characteristic of high voltage and 


low current. 


3. Electromagnetic Thrusters 


Electromagnetic thrusters (also known as magnetoplasma-dynamic or 


MPD thrusters) differ from electrostatic thrusters in that these utilize 


propellant gas that has been heated to a plasma state. There are several 
variations on this type of thruster, each of which requires its own unique and 
complex analytical model. Most, however, use plasma as part of a current- 
carrying electrical circuit which interacts with a magnetic field to produce 
thrust. 

Among the most interesting of the electromagnetic thrusters (and the 
focus of much of this report) are those known as Stationary Plasma Thrusters 
(SPTs), the design of which results in an electric field which is largely axial 
and points in the direction of the propellant flow, serving in turn to 
accelerate the ions in the exhaust to velocities as high as 16,000 meters per 
second, for a specific impulse of more than 1600 seconds (Brophy, 1992). 


Figure 3 shows a schematic diagram of an SPT. 
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Figure 3. Stationary plasma thruster schematic diagram (After Brophy, 1992) 
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Stationary plasma thrusters are often referred to as Hall thrusters or 
Hall accelerators, after the closed circular electron drift they exhibit between 
the cathode and anode, which is analogous to the Hall effect. Work on Hall 
accelerators originated in the United States in 1962 (Figures 4 and 5). 
However, U.S. research in this area had effectively ceased by 1970. By contrast, 
the former Soviet Union (present-day Russia) have employed stationary 


plasma thrusters for space applications since 1964, and their successes led to a 


reemergence of U.S. interest in the technology by the end of the last decade. 





Figure 4. Plasma accelerator using Hall currents (From Cann, 1962) 
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Figure 5. Annular magnetic Hall current accelerator (From Cann, 1964) 
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Figure 6. Stationary plasma thruster cross-section (From Day, 1995) 


External coil Ay es 


In 1991, a team of U.S. electric propulsion specialists visited Russia to 
experimentally evaluate the performance of a 1.35 KW stationary plasma 
thruster, designated the SPT-100 (Figure 6). Their examination verified that 
the actual performance of the thruster was close to the claimed performance: 
a thrust level of 80 mN at a specific impulse of 1600 seconds and an efficiency 
of 5090. Since that time, several SPT-100 thrusters have been made available 


to the United States for continued research and development at NASA's 


Lewis Research Center and Jet Propulsion Laboratory. Additionally, Space 
Systems/ Loral is currently flight-qualifying SPT-100 thrusters for north-south 


stationkeeping and orbit-raising applications on their satellites. 


B. ION ENGINE POWER THROTTLING 


Xenon ion propulsion systems have been proposed for use aboard 
interplanetary missions in the future, and NASA's NSTAR program will test 
the concept. Since the electrical power available to the thruster from solar 
energy will decrease with the square of the distance from the sun, these 
systems will have to be capable of power throttling across a wide range of 
inputs, a requirement that has not existed in the stationkeeping roles for 
which ion engines have been previously used. Accordingly, there is 
considerable interest in methods of throttling the output power of ion 
thrusters. Gridded ion thrusters such as NASA's 30-centimeter NSTAR 
engine may be throttled by varying discharge chamber voltage and current. 
These and Hall thrusters as well may also be throttled by varying the flow rate 
of the xenon gas propellant. It is this latter method which will be considered 


in detail in the rest of this report. 
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II. OVERVIEW OF CURRENT PROGRAMS 


A. NASA-JET PROPULSION LABORATORY ION PROPULSION 
INITIATIVES 


1. NSTAR 30-Centimeter Ion Engine 


The NSTAR program was established in 1993 as a means to assess the 
baseline capabilities of ion engine technology and confirm its potential to 
serve as a primary space propulsion system (Curran, 1995). The program is 
managed by the Jet Propulsion Laboratory, with NASA's Lewis Research 
Center responsible for developing the 30-centimeter xenon ion engine which 
has become the program's focus. Hughes Corporation is the prime contractor 
for the majority of the flight hardware, with the notable exception of the 
propellant feed system, which has been developed by Moog, Incorporated. 

NSTAR's electrostatic engine will provide the main propulsion system 
for NASA's New Millennium series of interplanetary missions, with the first 
launch scheduled for 1998. The 30-centimeter engine is considered ideal for 
these missions for a number of reasons. First, the engine itself is somewhat 
modular, allowing a spacecraft to use a main propulsion bank consisting of 
six hexagonally-arranged engines operating as a single unit (Figure 7). As 
higher power ion thrusters are developed, this arrangement is expected to 
eventually result in a segmented ion engine capable of power levels as high 
as 100 kilowatts, and specific impulses approaching 10,000 seconds, literally an 
order of magnitude greater than anything yet flown (Brophy, 1992). Second, 


the NSTAR engine is capable of being throttled across a wide range of output 
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powers, from 0.5 to 2.5 kW at efficiencies of 55-63%. This capability will allow 
the engine to vary its operation according to changes in the spacecraft's 


distance from the sun. 





Figure 7. Segmented Ion Thruster Arrangement (End View) 


2: Stationary Plasma Thruster 


As mentioned previously, the type of electromagnetic ion engine 
known as the stationary plasma thruster (SPT) was primarily developed in 
the former Soviet Union; in fact, these engines are often referred to as 
Russian Hall thrusters (RHT's). Since 1993, development of a 100-centimeter 
stationary plasma thruster, the SPT-100, has taken place under a joint 
venture, International Space Technology, Inc. (ISTI), founded by Space 
Systems/ Loral, the Russian Experimental Design Bureau Fakel and the 


Research Institute of Applied Mechanics and Electrodynamics (RIAME). 
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Performance, plume and EMI/ RFI tests were performed on the SPT-100 at 
NASA's Lewis Research Center in 1993, and a 5,730-hour endurance test was 
completed at the Jet Propulsion Laboratory in November, 1994 (Garner, 1995). 
The SPT-100 thruster body and xenon flow controller (XFC) are manufactured 
by Fakel. The thruster's power processing unit (PPU) was developed and is 
manufactured by Space Systems/ Loral, while the SPT-100's propellant 
management assembly (PMA) is built by Moog, Inc. 

Stationary plasma thrusters produce a specific impulse of about 1600 
seconds at an efficiency of more than 50%, and have demonstrated lifetimes 
in excess of 6000 hours. These characteristics make them ideal for north-south 
stationkeeping roles on geostationary spacecraft. Like the NSTAR engine, 
they may also prove to be suitable as main propulsion engines for 
interplanetary missions. Although stationary plasma thrusters have an 
extensive flight history on Russian spacecraft, their flight-qualification on 
western spacecraft has been delayed because of the need to redesign their PPU 
using western components and design practices in order to reduce subsystem 
mass and increase reliability (Day, 1995). 

One particularly unique feature included as part of the xenon flow 
controller on Fakel-built stationary plasma thrusters is an integrated 
“thermothrottle” used for fine control of the thruster’s propellant flow rate. 
The thermothrottle is a capillary tube flow control device which exploits the 
viscosity variation of xenon gas with temperature (Day, 1995). Chapters III and 
IV of this report will explore the concept of thermally throttling propellant 


flow in much greater detail. 
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B. ION ENGINE LIMITATIONS 


lon engines require an external source of electrical power, which is 
used to ionize and accelerate, and then to neutralize, the propellant. Thrust is 
produced only by propellant acceleration; the ionization and neutralization 
processes represent losses amounting to more than one-third of the input 
power. Other power losses may result from incomplete propellant ionization 
or premature neutralization, exhaust dispersion, and component heating. 

The power conditioning subsystem for an ion engine must typically 
provide outputs at two or three different grid voltages, two different cathodes 
and their heaters, and two or more electromagnets, as well as actuating power 
for isolation valves (Sutton, 1992). In some cases, ion engine PPUs have 
required over 4000 discrete parts and twelve power supplies in order to 
operate the thruster (Hamley, 1995); these PPUs can weigh more than the 
thrusters themselves, and can be more difficult to flight-qualify. 

The relatively low thrust levels produced by ion engines (typically less 
than 0.1 N) require that the thrusters operate for long periods of time in order 
to produce desired velocity changes. For example, using an ion propulsion 
system to raise a satellite from low earth orbit to geostationary orbit could 
require continuous thruster operation for as long as 200 days. Ion engines 
selected for interplanetary missions will be required to operate continuously 
and reliably for years without maintenance, a capability they have not yet 
demonstrated (Sutton, 1992). 

The need for reliable and wide-ranging power throttling of ion engines, 
especially those used in interplanetary roles, imposes additional constraints 


on the thruster’s power processing unit, which must be able to vary several of 
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the outputs mentioned above across the thruster's entire operating range. 
Power throttling also imposes additional requirements on the engine's xenon 
flow controller and propellant management assemblies, since the most 
commonly employed method of power throttling ion engines is to vary the 


thruster's propellant flow rate. 


c THESIS OBJECTIVES 


The thermothrottle device employed on Fakel-built Russian Hall 
thrusters represents a simple and lightweight means of throttling thruster 
propellant flow, since the device itself consists of little more than a capillary 
tube through which an electric current flows. However, current RHTs use the 
thermothrottle for only fine control of propellant flow, and the devices as 
currently employed require sophisticated propellant management assemblies 
to maintain propellant feed pressure within narrow limits. It is believed that 
the thermothrottle concept may, in an expanded role, offer much greater 
throttling capabilities for both the stationary plasma thruster and the NSTAR 
ion engine, and so may allow for future reductions in the complexity of both 
the power processing unit and propellant management assembly. The goal of 
this report, then, is to assess the theoretical capabilities and limitations of 
thermally throttling xenon propellant flow across a wide range of both flow 
rates and input pressures, and to explore ways in which thermal throttling 


may be used to reduce constraints on the PPU and PMA. 


15 


| 
| 


m — — ` 





16 


II. THEORETICAL BACKGROUND 


A. XENON PROPELLANT 


1. Advantages and Disadvantages 


The most desirable ion engine propellant would have a large atomic 
mass (which increases thrust); would be inert, inexpensive and easy to store 
and ionize; and would exist in a gaseous state throughout the propulsion 
system's storage and operating temperature ranges. Early ion engines used 
cesium or mercury as propellants; however, these elements have largely been 
abandoned due to their toxicity. Xenon, the noble gas with the highest atomic 
mass, is currently the preferred propellant for ion engines, although several 
argon-propelled systems are also under development. Rare and relatively 
expensive, xenon exhibits nearly ideal ion propulsion characteristics in every 
other aspect. A monatomic gas at room temperature, xenon is inert and easily 
ionized, and can be stored for years at moderate pressures without 


degradation. 


2. Thermoproperties 


Xenon exhibits essentially constant values of specific heat, gas constant 
and specific heat ratio. The other primary thermoproperties of gaseous xenon 
(enthalpy, entropy function, internal energy, thermal conductivity, sonic 


speed and viscosity) are functions only of temperature. Most significant to the 
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concept of thermally throttling xenon flow is the variation of the element's 
viscosity with temperature, which can be accurately calculated according to 


the following polynomial fit (Biblarz, 1981): 


u= 1.8918 x 10°° + 6.0506 x 10 T-8.1794x 10 T^ -2.0510x 10 ^T^ 
+ 2.0165 x 10°'°T* — 9.5234 x 10° T^ - 1.7663 x 10 ^T^ (3.1) 


Equations for the specific heat, gas constant, specific heat ratio, 
enthalpy, entropy function, internal energy, thermal conductivity, sonic 
speed and viscosity of xenon are given in Appendix A, as is MATLAB code 
for calculating the thermoproperties of xenon over a range of temperatures. 


Those properties are tabulated in Table A-1. 
B. FLOW REGIMES 


As described by Gombosi in his recent work on gas kinetics (1994), 
viscous gas flow falls into three different regimes: collision-dominated 
(Poiseuille) flow, -collisionless (free molecular) flow, and intermediate 
(transition or slip) flow. The different regimes can be characterized by the 
related dimensionless parameters of Mach number (M), Knudsen number 
(Kn) and Reynolds number (Re). The Mach number (Equation 3.2) is defined 
as the ratio is the bulk flow speed of the gas (v) to its sonic speed (a). The 
Knudsen number (Equation 3.3) is defined as the ratio of the mean free path 
separating collisions between molecules (A) and a characteristic linear size of 
the problem (L). The Reynolds number (Equation 3.4), which is defined as the 
ratio of the bulk momentum flux (Lup) and the viscosity of the gas (u), can 


also be expressed (Equation 3.5) in terms of the Mach and Knudsen numbers. 
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(3.2) 
» A 
n = 
L (3.3) 
p Lvp 
A= 
H (3.4) 
=M M 
Re =3 AAN (3.5) 


The Knudsen number has been found to be a good indicator of the 
collisional flow regime of a given problem (Gombosi, 1994). When Kr«1 (i.e., 
the mean free path is very small compared to the characteristic length of the 
problem), molecules undergo a large number of collisions while traveling 
any distances comparable to the scale of the problem, and the flow is collision- 
dominated. The collision-dominated regime can be easily analyzed using the 
methods of classical compressible fluid dynamics. This analysis can be further 
simplified if the flow velocity is small (the Reynolds number is less than 
about 2000), in which case the flow of the gas is considered to be laminar and 
is primarily controlled by its viscosity. 

At the other extreme, when Kn»1, the effects of intermolecular 
collisions are negligible and the flow is said to be collisionless. In this regime, 
molecules will collide with the walls of their container, but their collisions 
with each other may be ignored. Mathematical methods for the analysis of 


free molecular flow have been extensively developed in recent years. 
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The intermediate regime, when Kn=1, is mathematically the most 
difficult to describe and has been the least studied of the three flow regimes. 
There are currently no well-developed standard techniques for the general 
analysis of a transition flow problem. 

In the case of a propellant feed system using xenon gas at a pressure of 
two atmospheres and a temperature of 300 K, the mean free path is on the 
order of 10-8 meters. With a characteristic problem length of about 3 
millimeters (the inside diameter of the tubing used in a typical system), a 
Knudsen number of about 10° results, indicating that the problem is certainly 
one involving collision-dominated flow. Furthermore, it will be shown later 
that at the mass flow rates of interest (1-6 mg/sec), the gas velocity is low 
enough that the flow is entirely laminar and is therefore dominated by 
viscosity and temperature effects. 

Appendix B provides equations and MATLAB code relating isentropic 
flow properties to Mach number. These isentropic flow parameters are 


tabulated in Table B-1. 


C. FLOW DYNAMICS 


1. Fanno Flow 


Fanno flow is the adiabatic flow with friction of a compressible fluid in 
a constant-area duct. [Note: There are numerous outstanding texts which 
cover all aspects of Fanno flow and, as is considered in the next section, 
Rayleigh flow. The author relied primarily on Zucrow (1976) and Shapiro 


(1953) throughout this discussion of flow dynamics.] Since it is flow in which 
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friction alone is the driving force affecting the thermodynamic states of the 
fluid, Fanno flow is also referred to as simple frictional flow. Fanno flow is 
related to the bulk properties of the fluid through a friction coefficient f; in 
the case of laminar flow, the D'Arcy friction factor is a function only of the 


Reynolds number: 


_ 64 
Re (3.6) 
Equation 3.6 is valid when the Reynolds number is less than about 
2000, so it is useful at this point to confirm that this is in fact the case. Using 
the continuity equation for mass flow rate and the fact that the characteristic 
length of propellant feed problems like the ones to be considered in this 
report is equal to the diameter of the feed tubing, Equation 3.3 can be 


rewritten as follows: 


H uh (3.7) 


Using the limiting values of the maximum mass flow rate of interest 
(6 mg/sec) and (from Appendix A) the minimum xenon viscosity (which will 
occur at the minimum temperature of interest, 275 K), Equation 3.7 yields a 
maximum possible Reynolds number of 116, well below 2000. Accordingly, 
the gas flows of interest are all laminar, and Equation 3.6 is valid as a means 
of determining the friction factor. 

A plot of enthalpy (h) as a function of entropy (s) for a constant value 
of mass flux (G=pv) for the case of adiabatic constant-area flow with friction is 


called a Fanno line; such a plot is given in Figure 8. The Fanno line can be 
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divided into upper (subsonic) and lower (supersonic) branches, which, in the 
case of increasing frictional effects (either by increasing f or duct length, or 
decreasing duct diameter), converge at the maximum entropy point at the 
right. At this point, the flow velocity becomes sonic (M-1), and the flow is 
said to be choked. In other words, if the initial fluid velocity is subsonic (as is 
the case with propellant feed systems), increasing the effects of friction 
accelerates the fluid flow to sonic speed as a limit. Continuing to increase the 
frictional effects once the sonic limit has been reached will have the effect of 


causing the flow to move along a new Fanno line at a lower flow rate. 


Subsonic branch 






Supersonic branch 











Figure 8. A Fanno line in the h-s plane. 
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Appendix C provides equations relating Fanno flow properties to Mach 
number, along with MATLAB code used to calculate fluid properties over a 
range of Mach numbers. These Fanno flow parameters are tabulated in Table 
C-1. Appendix D provides the same information for the special case of 


isothermal Fanno flow. 


2 Rayleigh Flow 


Ideal Rayleigh flow is the frictionless flow with heat transfer of a 
compressible fluid in a constant-area duct. Since in this case heat transfer 
alone is the driving force affecting the thermodynamic states of the fluid, 
Rayleigh flow is also occasionally referred to as simple diabatic flow. The 
assumption that Rayleigh flow is frictionless is particularly valid in cases 
where the heat transfer is accomplished in such a short length of duct that the 
effects of friction can be considered negligible, as is likely to be the case with a 
thermal throttling system. The effect of Rayleigh flow on a fluid's bulk 
properties is normally determined through the use of basic heat transfer 
equations to determine changes in the stagnation temperature of the fluid; 
gas dynamics relations are then used to determine changes in the other 
thermodynamic properties of the fluid. 

As with Fanno flow, it is instructive to plot enthalpy (h) as a function 
of entropy (s) for a constant value of mass flux (G=pv) for the case of simple 
diabatic flow, producing a Rayleigh line such as is shown in Figure 9. The 
Rayleigh line too can be divided into upper (subsonic) and lower (supersonic) 
branches, which, in the case of increasing stagnation temperature, converge at 


the maximum entropy point at the right. At this point, the flow velocity 


23 


becomes sonic (M=1), and the flow may be considered to be thermally choked. 
Similarly to the case with Fanno flow, continuing to add heat once the 
limiting case has been reached has the effect of moving the flow to a new 
Rayleigh line at a lower flow rate. 

Appendix E provides equations relating Rayleigh flow properties to 
Mach number, along with MATLAB code used to calculate fluid properties 
over a range of Mach numbers. These Rayleigh flow parameters are tabulated 


in Table E-1. 


Subsonic branch 






h 


Supersonic branch 


Figure 9. A Rayleigh line in the /i-s plane. 
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3. Combined Effects of Friction and Heat Transfer 


Fanno flow assumes that there is no (or negligible) heat transfer, while 
Rayleigh flow assumes the same regarding friction; it may seem at first that 
the two are mutually exclusive. In the design of a thermally-throttled 
propellant feed system, however, it is necessary to consider the case where 
both heat transfer and friction are present. Shapiro (1953) presents a set of 
influence coefficients which relate the combined effects of changing area, 


friction and temperature on the Mach number of a flowing fluid: 








| 
dM? 2 + M | "m 
Ke A [area change] 
am E Y M 
+ Z ^p [friction] 
IA" D 
e£) fepe el Ac), 


dT, 
—-— [heating/cooling] 
IM" (3.8) 


Equation 3.8 shows that in the case of a constant-area feed system, 
provided that the heat transfer takes place within a relatively short length of 
tubing (the actual thermal throttle), it is reasonable to consider the effects of 
the heat transfer separately from those of the friction. In other words, by 
maximizing dT within a minimum dx (with dA=0), we can assume that only 
heat transfer affects the fluid within the thermal throttle itself. In all other 
sections of the feed system, both dT and dA are assumed to be zero or 


negligible, and the only effect is from friction. Figure 10 illustrates the Fanno 


(1—2)-Rayleigh (23)-Fanno (34) transitions such a system would exhibit. 
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Figure 10. Combined Fanno and Rayleigh Effects 


IV. ANALYSIS 


A. TEMPERATURE EFFECTS ON XENON FLOW 
1. Viscosity 


Recall from Equation 3.1 that the viscosity of xenon as a function of 


temperature (in the range of 120-1500 K) is given by: 


11.8918 x 10 ^4- 6.0506 x 10 T- 8.1794 x 10 !! T 22.0510 x 10 ^ P? 
-2.0165x 10 ^7'—9.5234 x 10 T^ - 1.7663 x 10 7 T^ (3.1) 


Figure 11 shows a plot of xenon viscosity vs. temperature; note from 
the detailed inset figure that in the temperature range of interest (275—375 K), 


xenon viscosity is essentially linear, and can be estimated from: 





u- [2.14640.00714 (T-275)]x10. || (4.1) 


2 Reynolds Number 


Equation 3.7 gave the Reynolds number in terms of mass flow rate, 
viscosity and tubing diameter. Combining Equation 3.7 with Equation 4.1, and 
assuming a typical tubing diameter of 3 millimeters, the Reynolds number 


can be expressed in terms of mass flow rate and temperature alone: 


4x 10° m 
£— 3x 2.146 +0.00714 (T— 275) (4.2) 


A 


Temperature (K) 


1500 


1000 


500 


3 4 5 6 
Viscosity (N-sec/m?) 


Figure 11. Xenon viscosity as a function of temperature. 
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Figure 12 shows the Reynolds number as a function of temperature 
determined using Equation 4.2 for several mass flow rates in the range of 


interest for propellant feed systems. 
3. D'Arcy Friction Coefficient 


Equation 3.6, which gave the D'Arcy friction factor as a function of the 
Reynolds number for laminar flow, can be combined with Equation 4.2 to 


relate the friction factor to temperature and mass flow rate: 


487 2.146+0.00714(T-275) 
10° m (4.3) 


Figure 13 presents the D'Arcy friction factor as a function of xenon 
temperature for several mass flow rates of interest. (Note: several texts use 


the Fanning friction coefficient, which is equal to f/4, instead.) 
4. Mass Flow Rate 


Temperature effects on the flow of xenon occur as a result of changes in 
the friction factor. The effect of friction varies directly with the length and 
inversely with the diameter of the duct. Tables of Fanno flow parameters 
(such as Table C-1) account for these variations by including the length-to- 
diameter ratio in their calculations, which result in a critical length parameter 
AfL / D. The effect of friction (and therefore temperature) can be seen by 


comparing isentropic choking flow to frictional choking flow, as a function of 
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4fL/ D. Figure 14 presents such a comparison for y-1.67. The figure is derived 


from the fact that, for a choked flow, 


IP:A* 


p — op ue KC 
7 + ] 
WRT, d (4.4) 








(Note: wherever used in this report, the superscript * refers to the critical, or 
sonic, condition; the subscript t refers to the stagnation state.) Since, for both 
Fanno and isentropic flow, the stagnation temperature is constant, the ratio of 
the Fanno and isentropic flows can be expressed simply as the ratio of the 
stagnation pressures, which is one of the Fanno flow parameters given in 


Appendices C and D. 


B. BASIC SYSTEM ANALYSIS 


It should be apparent from the preceding section that, in the case of a 
thermal throttling system, mass flow rate and frictional losses are mutually- 
dependent, a fact which complicates the analysis of such systems. Further, 
traditional simplified gas dynamics methods, which generally assume some 
knowledge of conditions at each end of a duct (or a given mass flow rate 
through it), are insufficient for the analysis of a system in which the fluid 
undergoes a series of transitions between the points of known conditions, 
such as is shown in Figure 15. Analysis of these systems is likely to require an 
iterative process based on some simplifying assumptions. Such a process is 


presented next. 
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Figure 15. Basic thermothrottle system. 


1. General Method 


Consider the basic thermothrottle system shown in Figure 15, in which 
a large tank of xenon at constant pressure and temperature feeds a long 
constant-area duct with a short thermal throttle between the ends (Stations 
2—3). The duct is assumed to be well-insulated, and exhausts to a vacuum. 
The duct inlet (Station 1) is considered to behave as an isentropic nozzle. The 
temperature at the outlet of the thermal throttle (Station 3) is known. The 
flow through the duct is considered to be isothermal Fanno (Appendix D) 
from Stations 12 and 34, and Rayleigh (Appendix E) from Stations 23. 

The general method to be used is to determine the maximum 
isentropic flow rate through the duct, and to use this flow rate as a means of 
estimating Mach number and friction coefficient at Station 1. Fanno flow 
analysis is then used to estimate conditions at Station 2, and Rayleigh flow 


analysis is used to estimate those at Station 3. The final section of the duct is 
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considered to be choked at Station 4, and the conditions at Station 3 are used 
to determine a more accurate friction factor. The problem is then worked back 
to Station 1 to establish more accurate Mach numbers at each station. These 
values are used to determine final values for the pressure, density and flow 
rate at Station 3. The problem can then be reworked for additional 
combinations of inlet pressure and throttle temperature, or for variations in 
duct length and/or diameter. 

Equations 4.5 through 4.8 allow determination of the maximum 


isentropic flow rate as a function of tank pressure and temperature. 





Py 
RT 
(4.5) 
2 

y*- um RT, 
Ye : (4.6) 

1 

E E 

mor 
EE (4.7) 
Misen= P A'V" (4.8) 


Once an inlet flow rate has been determined, Equation 4.3 can be used 
to calculate the friction coefficient in the duct. Equation 4.9 can then be used 


to determine the Mach number at any station: 








D WM (4.9) 


Mach number determination is normally accomplished through the 
use of tabulated values of the critical length ratio 4fL*/ D as a function of inlet 
Mach number. In this case, however, it is desired to allow the direct 
calculation of the Mach number for a given value of the critical length ratio. 
Equation 4.9 cannot be solved explicitly for M, but a series expansion may be 


employed in place of the logarithmic term (Equation 4.10): 
l 3 


4fL* _ 1— yM” ; i 
PM —— — 
D m 2 (4.10) 








The Mach number is considered to be low enough at Stations 1—3 that 


the higher-order terms in Eguation 4.10 can be neglected, resulting in: 


a. La 
D` yM’ 








va 1)=| Te 





1 
YM” (4.11) 
It has been determined empirically that the slight modification of 

Equation 4.12 provides a more accurate curve fit at moderate Mach numbers 
(below 0.5). Figure 16 shows the correlation between the exact function 
(Equation 4.9) and the estimate from Equation 4.12. Equation 4.13 gives the 
desired solution of Equation 4.12 for the Mach number as a function of the 


critical length ratio: 














4fL* [0.9 
E -| 2 ryM^-2 
(4.12) 
— y 425y^100y 4-10 _ 4fL 
M= Wiad 10 where y = D (4.13) 


36 


4500 ——— Actual (non-linear) function 
— — — Linear approximation 


4000 


3500 


2 


o 
2 
a 
a 
£ 
D 
E 
o 
= 
S 
=D 
E 
o 


Critical Length Ratio 
on 
e 
o 


Š 


1500 


1000 0.25 
Mach Number 





O 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
Mach Number 


Figure 16. Critical Length Ratio linear approximation. 


The Mach number at Station 1 can now be used (in Equations 4.14 
through 4.16) to determine the density at Station 1, along with a revised 


frictional flow rate: 


ya p 
=p,|1+——M 
A nl EEN (4.14) 
p - p, JY M, (4.15) 
M pene RAEV" (4.16) 
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The solution process is then continued using standard isothermal 


Fanno flow techniques from Stations 12 and Rayleigh flow techniques from 


Stations 23. The result will be revised values for the friction coefficient in 
both Fanno flow sections, which will then allow for determination of final 
values for the Mach number at each station, and a final value for the mass 
flow rate through the duct. Appendix F contains MATLAB code (flowcalc.m) 
written to automate this process, and to allow for variations in the inlet 


pressure and in the length and/or diameter of the duct. 


2: A Numerical Example 


To illustrate the process, consider a system such as the one shown in 


Figure 15, with initial conditions as given in Table 1 below: 





275 K 
2 







Table 1. Sample system initial conditions. 


The MATLAB program flowcalc.m (Appendix F) is used to analyze the 
system. The code determines a maximum isentropic flow rate of 
389 mg/sec, and an unthrottled (275 K throughout the duct) Fanno flow rate 


of 129 mg/ sec, with an initial critical length ratio (4fL*/ D) of 16.6. Note that 
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this value, when used as the input to Figure 14, suggests a flow ratio of about 
0.3; the actual calculated ratio is 0.33. The throttled flow rate was calculated to 
be 5.34 mg/sec, with a critical length ratio of 41 in the first Fanno section and 
970 in the second. These values suggest a flow ratio of about 
(0.25)(0.05)=0.0125; the actual calculated ratio is 0.014. The analysis was then 
conducted over the entire temperature range of the thermothrottle (275-375 
K) for several different pressures. Figure 17 shows a plot of the resulting 


calculated mass flow rates as a function of temperature. 


(e COMPARISON WITH REAL SYSTEMS 


There is very little experimental data available on any of the thermal 
throttling systems currently in use. Most of these systems are based on 
former-Soviet designs which have been used successfully for many years, but 
which require very tightly controlled inlet conditions and exhibit very 
narrow operating ranges. Current ion engine systems using thermal throttles 
have been designed around these existing throttles. The author is unaware of 
any case in which a thermal throttle has been operated independently of its 
engine for the purpose of developing a full set of operating range data—in 
other words, no one has yet opened up the "black box" to see exactly what is 
inside. As a result, it is difficult to perform any quantitative comparisons of 
the predicted flow rates obtained using the methods documented in this 
report with actual flow rates observed in comparable real-world systems. 

The author was given limited access to two Russian ion engine 
propellant feed system designs, although in neither case was sufficient 


information available. Significantly, however, the shapes of the flow rate 
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Figure 17. Throttled flow as a function of temperature 


for several input pressures 


curves in Figure 17 are well-supported by the available data, suggesting that 
the approach used in this report is valid. Other information gained from 
studying real-world systems is that they are fairly compact («1 meter from 


pressure regulator to engine) and that they operate at inlet pressures of 2-3 
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atmospheres. These systems can be modeled using the MATLAB routine 
flowcalc.m by varying the input pressure and length parameters, resulting in 
reasonable agreement with existing systems' measured flow rates. 

It should be apparent that the combination of increased inlet pressure 
and reduced duct length will result in higher flow rates; flowcalc.m can easily 
verify that fact. Existing systems compensate for their requirements of 
compactness and higher operating pressures through the use of flow 
restrictors. There are currently two types of flow restrictors in use in ion 
engine propellant feed applications. The first, generally referred to as a flow 
resistor, is shown in Figure 18. Flow resistors are simply sections of duct with 
reduced areas; they affect flow by increasing the L/ D term in the critical length 
ratio, and can be modeled linearly for laminar flow using a similarity rule 


developed by Lee (1974). 


Dn 


| m 
EE . 


Figure 18. Flow resistor. 


The second type of flow restrictor consists of a porous metal plug, 
occasionally referred to as a sintered restrictor. These are generally designed 
for a specific combination of gas, input pressure and flow rate. Preliminary 
work by Martin (1974) using nitrogen gas indicates that, once sufficient data 
are obtained for a given plug and system combination, it becomes relatively 


easy to predict the performance of the system under a wider variety of 
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conditions provided the pressures at the inlet and outlet of the plug are 
known. That requirement presents a considerable challenge, however, to an 
analysis based on Fanno and Rayleigh techniques, since a porous plug may be 
considered as an isenthalpic device which cannot be simply modeled as an 
equivalent L/ D. Figure 19 shows the enthalpy-entropy diagram which may 
result from the addition of a porous plug at Station 4 in the basic 
thermothrottle system of Figure 15. Modeling such a system will likely 
require some knowledge of the temperature and pressure conditions at 


Station 4, and so will require some degree of instrumentation. 
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Figure 19. Combined process with isenthalpic flow restrictor. 


Another area of interest in the design of current and future thermally- 
throttled ion engine propulsion systems is the possibility of using the thermal 


characteristics of the propellant as a substitute for expensive and failure-prone 
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pressure regulators. Manipulation of the various input parameters in 
flowcalc.m suggests that a thermothrottle system can be designed to operate 
over a wider range of input pressures than that used in current systems. 
However, complete elimination of the need for a pressure regulator would 
require a system capable of operating under input pressures varying from 
1-100 atmospheres. While this could be accomplished using a bank of isolable 
flow restrictors (which would allow for an extremely large change in L/D 
over the life of the system), the solenoid valves and tubing connectors 
required in such a design would add both mass and complexity to the system, 
and so this likely not an acceptable alternative to the commercially-available 
pressure regulators currently in use. 

Real-world thermal throttling systems generally regulate flow rate as a 
function of the electrical current supplied to the throttle. This current is used 
to heat the gas in one of two ways, either through use of an external heating 
element or by direct current flow through the walls of the throttle; in both 
cases, heat is produced through IFR losses. Current flow can be related to gas 


temperature through a basic heat transfer relationship: 


Nu pk 


= P R =|- —— | (T. -T 
| D | Ty) (4.17) 


where Nu, is the Nusselt number, k is the thermal conductivity of the gas, D 
is the tubing diameter, and Tọ and Tp are the wall and bulk temperatures, 
respectively. For the uniform heat flux, laminar flow case of a thermal 
throttle, Nu, is a constant (4.364), and Tp (which equates to the stagnation 
temperature of the gas) varies linearly with q (and so with 12). As a first 


approximation, it seems that gas flow (which, from Figure 17, varies inversely 
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with temperature) will, in a lossless system, vary inversely with the square of 
the current. However, in the uninsulated systems the author has observed, 
gas flow has tended instead to vary as the inverse of the current. 

It will be difficult to predict the exact relationship between throttle 
current and gas flow rate without knowing which heating method is 
employed, as well as how effectively the system is insulated, since those 
factors affect how much of the supplied current is actually used to heat the 
gas. Also, it must be realized that ground testing is subject to free convection, 
which depends on the orientation of the heated tubing with respect to gravity; 
this effect further complicates determination of the current and flow rate 
relationship. Accordingly, this report has focused on the more well-defined 


correlation between gas temperature and flow rate. 


V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


This report has presented a possible method for predicting the 
performance of thermally-throttled ion engine propellant feed systems, based 
on an adaptation of the traditional gas dynamics approach. While this 
method appears to agree qualitatively with observed system behavior, it is 
incomplete in that the means to deal with system components such as porous 
plugs have not yet been developed. 

The results obtained in the previous chapter suggest that a single 
thermothrottle system, while capable of operating across a wider range of 
inputs and outputs than is currently seen, is nonetheless unlikely to allow for 
the complete replacement of the pressure regulators now in use. However, 
the thermothrottle system does allow for limited variation in input pressure 
while still delivering the desired flow rates, as opposed to the rigid control of 
input pressure required in current systems. This flexibility may at least allow 
for some simplification or reduction in demands on the systems' power 


processing unit, without much of an additional mass or power requirement. 


B. RECOMMENDATIONS 


Future study of the thermothrottle concept should focus on an 
empirical analysis of an operating system as a means of validating or 


contradicting the approach presented in this report. A simple basic system 
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such as that depicted in Figure 15 could be easily and inexpensively 
constructed, and would require only the installation of appropriate 
instrumentation to assess the accuracy of this report's methods. Such 
instrumentation should include the means to directly measure gas flow rate, 
such as through an instrumented orifice or hot wire flow sensor, along with 
pressure and temperature sensors at the inlet and outlet of the thermothrottle 
and any system flow restrictors. Given the expense and relative scarcity of 
xenon gas, argon is suggested as an appropriate substitute, since its 
characteristics and properties are similar to those of xenon. In the MATLAB 
program flowcalc.m, for example, only the gas constant R would need to be 
changed to allow for the analysis of argon in place of xenon. 

Alternatively, instrumentation could be installed on one of the 
thermothrottle systems currently on loan to NASA (which already have the 
means to measure flow), and direct measurements of the gas temperatures at 
the throttle inlet and exit, as well as the pressures across any flow restrictors 
in the system, could be integrated into ongoing ion engine test programs. 
Since knowledge of the temperature rise across the throttle (along with 
knowledge and/or control of heat losses in the system) allows determination 
of the relationship between current and temperature, and knowledge of the 
pressures at the inlet and exit of flow restrictors allows the remainder of the 
system to be analyzed according to the Fanno and Rayleigh relationships, 
these measurements would provide the data necessary to complete the 


thermothrottle performance prediction model. 


APPENDIX A. XENON THERMOPROPERTIES 











EQUATIONS 
1. Specific heat capacity (constant pressure) 
— 158 "A REN 
c,= 158. CS 
2. Gas constant 
R J 
R EM 63.323 kee K 
3. Specific heat capacity (constant volume) 
=c,—R=95.093 : 
E 
4. Specific heat ratio 
7 1.67 
E&E C, ed, 
5. Enthalpy 
kJ 
Gel c,dT—c,T E 
6. Internal energy 
kJ 
uH. ed E 
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7. Entropy function 


dT kJ 
me = catar [| 


8. Acoustic speed 
m [*] 


9: Viscosity 





[= 1.8918 x 107° + 6.0506 x 10°* T+ 8.1794 x 10 T^-2.0510x 10 ^T | Nes 
t2.0165:x 107^ T*—9 5234 0'r “O | m’ | 


Si Thermal conductivity 





y W 
k 21.3554 x 10 2.0398 x 10^ T—5.5381 x 10 ^ T? | | 


meK 


MATLAB CODE 


thermoproperties.m 
Thermal Properties of Xenon vs. Temperature 


Variables (all in SI units): 


= Temperature 

- Specific Heat at Constant Pressure 
- Specific Heat at Constant Volume 
- Gamma (Specific Heat Ratio) 

- Gas Constant 

Enthalpy 

- Internal Energy 

- Entropy Function 

- Viscosity 

- Thermal Conductivity 

- Local Speed of Sound 


e c oe 
| 
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ae 120 3110: 100016 
T2 = 1050550: 1500] ; 
EL T T2]; 

mo 298.15; 


X = ones(size(T)); 


CP - 158.416*X/1000; 


"m. 67*X; 

ies .51434/131.29; 
“U—CP- R; 

fe CP.*(T); 

pa ho (RT); 


PHI — CP.*(log(T)); 


ni  E-c) (1189178728 + (6.0506328E-2)*T + (8.1793523E-5)*(T.%*2) 
= (2.051E-7)*(T.^3) + (2.0164963E-10)*(T.^4) - (9.5234E-14) * (T.^5) 
1:8 76626569E-17) ^(T.^6)); 


5E es D4 26R-6 + (2.03984915E-5)*T - (5.53807454E-9) * (T.^2) ; 


EEssuart (g*1000*R- *T); 


TEE Du: BI mk; al; 
fid - fopen('thermoproperties.txt','w'); 
fprintf (fid, '%2.2f\t %4.4f\t $4.4fNC $4. 4fNC $4.4eNC $4.4eNC ... 


Si An! Y); 
fclose(fid); 
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Table A.1. Thermodynamic Properties of Xenon in SI Units 


T [K] 


120: 
1302 
140. 
T505 
160. 
170 
T80: 
T907 
200. 
ZU 
220 
DU 
240. 
230. 
260. 
270. 
280. 
2907 
3007 
SiG: 
3207 
290; 
340. 
3905 
3607 
SIE 
380. 
390. 
400. 
410. 
420. 
430. 
440. 
450. 
460. 
470. 
480. 
490. 
500. 
S10. 
320. 
930; 
540. 
Su. 
90D. 
559. 
nou 
SUD. 
600. 
610.1 
620. 


kJ | 

kg u 
20099 
.9941 
1782 
. 7624 
. 3466 
023507 
.5149 
, 0290 
25952 
.2674 
ENE 
2429 
70198 
. 6040 
. 1882 
LE 


kJ kJ 

kg | ^|kg*K 
oa 0.7584 
250151 0.7711 
a 0.7828 
2032 0.7938 
STE 0.8040 
. 1650 0.8136 
Ss 0.8226 
-0667 0.8312 
. 0176 @ 
. 9685 (OT 
. 9193 0.8544 
. 8702 0:8615 
821 0.8682 
3120 058727 
. 1229 0.8809 
265737 0.8869 
26246 0.8926 
ub 55 0.8982 
5264 079036 
sd 0.9088 
. 4281 0.9138 
3790 975 
.3299 0.9234 
2808 0.9280 
<2317 0.9325 
moss 0.9368 
mp4 0.9410 
-0843 0.9451 
20552 0.9491 
aol 0.9531 
. 9369 OSOS 
.8878 0.9606 
-8387 0.9642 
. 7896 0.9678 
dos 0.9713 
. 6913 0.9747 
26222 0.9780 
5931 0.9813 
.5440 0.9845 
. 4949 0.9876 
-4457 0.9907 
E5956 0:992 
23475 0.9967 
. 2984 0. 9995 
. 2493 1.0024 
2001 1.0053 
ISO 1-0080 
-TOTS BROTO 
20528 qa 
-0036 1.0160 
-9545 1.0186 
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| H m? 
.0015e-05 
.0743e- 05 
.1476e-05 
.2211e-05 
.2949e-05 
.3689e-05 
.4431e-05 
.5174e-05 
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.7403e-05 
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.9625e-05 
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.4004e-05 
.4723e-05 
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. 6147e-05 
. 6853e-05 
.7555e-05 
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.8944e-05 
.9632e-05 
.0314e-05 
.0992e-05 
.1664e-05 
.2332e-05 
.2994e-05 
.3651e-05 
.4303e-05 
.4950e-05 
.5591e-05 
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.7485e-05 
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.8721e-05 
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. 0538e-05 
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dd 


.3694e-03 
.5596e-03 
. 7486e-03 
.9365e-03 
.1233e-03 
. 3090e-03 
.4937e-03 
.6771e-03 
. 8595e-03 
. 0408e-03 
.2210e-03 
. 4000e-03 
.5780e-03 
. 7548e-03 
.9306e-03 
.1052e-03 
.2787e-03 
.4512e-03 
.6225e-03 
.7927e-03 
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.9695e-03 
.0109e-02 
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6803 
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0817 
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. QUEE 
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149. 
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4358 
0273 
5343 
9625 
3163 
6022 
8223 
98 
0819 
129 
121 
0661 


8157 
625 
3955 
1226 
8141 
4694 
0900 
6/3 
2324 
7563 


. 2 
2158 
2188 
2208 
2228 
2250 
2208 
2238 
. 2428 
234% 
2568 
2388 
241. 
243. 
245. 
247. 
249. 
23] 
258. 
256. 


7161 
1599 
5643 
9488 
3082 
6430 
9542 


5082 
1828 
9754 
1780 
3606 
9299 
6682 
7942 
9022 
9929 
0661 


Table A.1. Thermodynamic Properties of Xenon in SI Units (cont.) 
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3862 
9704 
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2212 
4754 
0595 
6437 
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8120 
3962 
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1486 
2328 
5: 10970, 
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6536 
25 16 
8219 
4061 
3902 
5744 
1586 
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2576 
1784 
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0200 
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8616 
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7037 
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kJ 
u kg 
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T380 
1/089 
-6598 
OLOT 
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. 4633 
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.2668 
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DE 
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0217 
29723 
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EE 
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6774 
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235000 
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.4318 
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139536 
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11262 
ES Jal: 
20579 
. 8423 
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SE 
ASS 
203999 
. 6143 
29557 
2182 31 
TD 
26319 


| ol 


PREP PPP PPP BP PEP PPP PPP B. B. S. BS. p. S. B. B. BB. p. B. S. B. B3. BS. B3. p. ES. ES. ES. B. p. p. p. p. p. [B HB HS 





kg*K 


20211 
20256 
.0261 
20255 
20509 
505.32 
EUIS S5 
-0378 
.0400 
. 0423 
.0444 
.0466 
.0487 
-0508 
20529 
.0549 
29570 
10569 
20609 
Es. 
. 0648 
. 0667 
. 0686 
. 0704 
20022 
.0740 
320958 
076 
. 0794 
S08 TL 
. 0828 
. 0845 
20862 
0578 
20895 
OS LI 
20927 
.0943 
02.0 
21094 
.1164 
202.32 
96 
2359 
.1418 
.1476 
11532 
iS 


kJ 
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o o - sl sl sl sl CDN OY OY OY OY OY OY OY OY OY U1 U1 U1 U1 U1 ua AMAN U1 U1 Ul Ul Ul Ul Ul U! UI UI Ul Ul Ul UI GGG VV eoe 





Nes 

| H m 

.4609e-05 
.5172e-05 
.5731e-05 
.6286e-05 
.6836e-05 
.7383e-05 
.7925e-05 
.8464e-05 
. 8999e-05 
.9530e-05 
.0057e-05 
.0581e-05 
.1101e-05 
.1618e-05 
.2132e-05 
.2642e-05 
.3149e-05 
.3653e-05 
.4154e-05 
.4652e-05 
.5147e-05 
.5639e-05 
. 6129e-05 
. 6616e-05 
.7100e-05 
.7581e-05 
.8060e-05 
.8537e-05 
.9011e-05 
.9482e-05 
.9951e-05 
. 0418e-05 
. 0882e-05 
.1344e-05 
.1804e-05 
.2262e-05 
.2717e-05 
.3170e-05 
54036-05 
.7581e-05 
.9706e-05 
.1778e-05 
.3796e-05 
.5765e-05 
.7689e-05 
.9579e-05 
.1452e-05 
.3331e-05 


a 


L3 H P H P H P P H HP H H H H H BH HB HP H H H H H PH H pP| H H H pP| H BH HB HP H BI pmpmpbmppmpmpegmpermppnmp 





Ee? 


.0654e-02 
.0788e-02 
.0921e-02 
.1052e-02 
.1182e-02 
.1312e-02 
.1440e-02 
.1567e-02 
.1693e-02 
.1817e-02 
.1941e-02 
.2064e-02 
.2185e-02 
.2305e-02 
.2425e-02 
.2543e-02 
.2660e-02 
.2776e-02 
.2891e-02 
. 3004e-02 
.3117e-02 
.3228e-02 
.3339e-02 
. 3448e-02 
. 3556e- 02 
. 3663e-02 
. 3769e-02 
. 3874e-02 
. 3978e- 02 
.4081e-02 
.4182e-02 
.4282e-02 
.4382e-02 
.4480e- 02 
.A5/7e-02 
.4673e-02 
.4768e-02 
.4862e-02 
.5314e-02 
.5739e-02 
.6136e-02 
. 6505e-02 
. 6846e-02 
.7160e-02 
.7446e-02 
. 7705e-02 
. 7/935e-02 
.8138e-02 


259: 
260r 
252 
264. 
206: 
268. 
DUE 
DYD 
274 
Du 
2E 
TO 
2 
20a. 
Z GS 
Zu 
289. 
290. 
292 
294. 
296: 
28. 
2293 
a Que 
303 
3051 
306. 


S10: 
Sq 
DES 
315 
Bos 
Soe 
520% 
DNE 
DU gd 
GE 
EE 
34]. 
348. 
356 
DOS. 
370. 
Sr 
384. 
ao. 
395. 


1229 
1635 


dL 


1972 
dS HEP, 
1703 
i350 


20853 
20220 


9450 
8547 
75 is 
9552 
5066 
3656 
2127 
0479 
8716 
5539 
4850 
2192 
0547 


.8236 


5821 


. 3304 


0687 
7971 


NI 


22951 
9250 


25157 
32073 


9700 
EE 
2991 
9358 
5742 
2043 
239 
(07870 
7428 


22435 


3595 
7900 
8335 
7869 
5979 
2923 





APPENDIX B. 


ISENTROPIC FLOW PROPERTIES 








EQUATIONS 
1. Dimensionless velocity 
[I ew 
Y 2+(y-1)M” 
D: Temperature ratio 
T YM i 
DR 14M | 
3. Pressure ratio 
E Vo T 
pu 1+ > M | 
4. Density ratio 
-1 
pr mo Mir 
p. 2 
5. Area ratio 


> 





> 


3o 


B. MATLAB CODE 


Isentropic.m 
Isentropic flow parameters 
Variables: 

M - Mach Number 


k - Gamma (Specific Heat Ratio) 
MSTAR - Dimensionless Velocity 


GG 00 o? de o9 o? oe o? o9 oe o9 oe oe oe oe 


TRAT - Stagnation Temperature Ratio 
PRAT - Stagnation Pressure Ratio 
RHOTRAT - Stagnation Density Ratio 
ARAT - Critical Area Ratio 
= O ON E 


ic xg og 
MSTAR - sqrt(((k-1).*M.^2) ./ (2* (k-1) . *M.^2)) ; 


TRAT 


dey CE EICH E 


PRAT 


(L-((k-1)72). *M. ^2) ea 

REORAT - (del ka ee ENDE 

ARAT = (1./M) .*((2/(k+l1)) . * (14+ ((k-1) 2). MEE NEE: AE 
Y- [M; MSTAR ; TRAT; PRAT; RHORAT; ARAT] ; 

fid = fðpen('isentropic. tt we 


fprintf (fid, '%2.2£\t %4.4e\t %4.4e\t %4.4e\t $4.4eyt $4.4e'n',Y); 
fclose(fid); 


54 


Table B.1. Isentropic Flow (y=1.67) 


cococccccococc. cc oc Ooo 05 020000050 GC: oc 02 (0602. 65 (60 6500 0 2OA2X: Oc CS DS cc 0 C9 C2 CO CO C9) oOo OQO 


un au u ui ui UI i> i> aS aS aS aS aS aS UU. UU. UU) UJ. UJ. CJ. CJ GJ GJ NO NO NO NO BO BO BO NO NJ) IP. ED. ID. D. Bj. ES. p. OB. BS. X0 O0 OY Ui iS UJ NO. HP. o 


M* 


.0000e+00 
.1554e-02 
.3107e-02 
.4657e-02 
.6204e-02 
.7747e-02 
. 9284e-02 
. 0813e-02 
.2335e-02 
. 0385e-01 
.1535e-01 
.2684e-01 
.3832e-01 
.4978e-01 
.6123e-01 
.7266e-01 
.8408e-01 
.9548e-01 
.0686e-01 
.1821e-01 
.2955e-01 
.4087e-01 
.5216e-01 
.6342e-01 
.7466e-01 
.8588e-01 
.9706e-01 
.0822e-01 
:1935e 0T 
.3045e-01 
.4152e-01 
.5255e-01 
.6355e-01 
.7452e-01 
.8545e-01 
.9635e-01 
. 0721e-01 
.1803e-01 
.2881e-01 
.3956e-01 
.5026e-01 
.6092e-01 
.7154e-01 
.8212e-01 
.9266e-01 
. 0315e-01 
. 1360e-01 
.2401e-01 
.3436e-01 
.4468e-01 
,5494e-01 


LO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO MO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO AO XO XO XO XO XO XO XO AO MO AO AO AO VO VO P 


7 


T 


t 


.0000e+00 
.9997e-01 
.9987e-01 
.9970e-01 
.9946e-01 
.9916e-01 
.9880e-01 
.9836e-01 
.9786e-01 
.9729e-01 
.9666e-01 
.95966-01 
.9520e-01 
.9437e-01 
.9348e-01 
.9252e-01 
.9150e-01 
.9041e-01 
.8926e-01 
.8805e-01 
.8678e-01 
.8544e-01 
.8404e-01 
.8259e-01 
.8107e-01 
. 7/949e-01 
.7786e-01 
.7616e-01 
.7441e-01 
.7260e-01 
.7073e-01 
.6881e-01 
.6683e-01 
.6480e-01 
.6272e-01 
.6058e-01 
.5839e-01 
25615e-01 
.5386e-01 
.5152e-01 
.4913e-01 
.4669e-01 
.4420e-01 
.4167e-01 
.3909e-01 
.3647e-01 
.3381e-01 
.3110e-01 
.2835e-01 
.2555e-01 
.2272e-01 


CO CO CO 00 OO 00 00 00 00 OO OO OO OO OD CO XO XO XO XO XO XO XO XO XO XO XO XO XO NO XO XO XO XO XO XO XO MO XO XO XO VO XO XO XO VO LO VO VO MO MO ES 


Jg 
Jg 


Le 


P 


t 


.0000e+00 
.9992e-01 
.9967e-01 
.9925e-01 
.9867e-01 
.9792e-01 
. 9700e-01 
.9592e-01 
. 9468e-01 
.9327e-01 
.9170e-01 
.8997e-01 
.8808e-01 
.8603e-01 
.8382e-01 
.8146e-01 
.7894e-01 
.7627e-01 
.7345e-01 
.7048e-01 
.6737e-01 
.6411e-01 
.6070e-01 
.5716e-01 
.5348e-01 
.4966e-01 
.4571e-01 
.4163e-01 
.3742e-01 
.3309e-01 
.2864e-01 
.2406e-01 
.1937e-01 
.1456e-01 
.0964e-01 
.0462e-01 
.9949e-01 
.9425e-01 
.8892e-01 
.8349e-01 
.7797e-01 
.7236e-01 
.6666e-01 
.6088e-01 
. 5502e-01 
.4908e-01 
.4307e-01 
.3699e-01 
.3084e-01 
.2462e-01 
.1835e-01 


CO OO OO CO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO MO XQ XO WO XO XO VO VO VO H 


Us 
p, 


.0000e+00 
.9995e-01 
.9980e-01 
.9955e-01 
.9920e-01 
.9875e-01 
.9820e-01 
.9756e-01 
.9681e-01 
.9596e-01 
.9502e-01 
.9398e-01 
.9284e-01 
.9161e-01 
.9028e-01 
.8885e-01 
.8734e-01 
.8572e-01 
.8402e-01 
.8222e-01 
.8033e-01 
.7835e-01 
.7628e-01 
.7412e-01 
.7188e-01 
.6955e-01 
.6713e-01 
.6463e-01 
.6204e-01 
.5938e-01 
.5663e-01 
.5381e-01 
.5090e- 01 
.4792e-01 
.4487e-01 
.4174e-01 
.3854e-01 
.3526e-01 
.3192e-01 
.2851e-01 
.2503e-01 
.2149e-01 
.1788e-01 
.1421e-01 
.1047e-01 
.0668e-01 
.0283e-01 
.9893e-01 
.9497e-01 
.9095e-01 
.8688e-01 


Lp) p ppprpppmp»mpm»pmp»fpmP.?p»pmp pup pm p pt) NN N NN NN N IN IN NN) U U UJ C GU iis iS iS UI UI OY - OO X0 ID HD ES INO UT 


A 
A 





oo 


.6235e+01 
.8123e+01 
.8/55e+01 
.4073e+01 
.1265e+01 
.3943e+00 
.0593e+00 
.0589e+00 
.2817e+00 
.6607e+00 
.1533e+00 
.7310e+00 
.3744e+00 
.0692e+00 
.8052e+00 
.5747e+00 
.3718e+00 
.1919e+00 
0313e 100 
.8871e+00 
.7571e+00 
.6392e+00 
.5319e+00 
.4339e«00 
.3440e+00 
.2614e+00 
.1852e+00 
.1147e+00 
.0494e+00 
.9886e+00 
.9321e+00 
.8794e«00 
.8301e+00 
.7839e«00 
.7406e«00 
. 7/000e00 
.6618e+00 
.6258e+00 
.5919e+00 
.5599e+00 
.5297e+00 
.5011e+00 
.4740e«00 
.4484e400 
.4242e400 
.4011e-«00 
.3793e«00 
.3585e+00 
.3388e+00 
.3201e+00 


Table B.1. Isentropic Flow (y-1.67) (cont.) 


Soom OO OOO COCO COO OO O10 CO CeO. Ooo. O 053 03 0 (CO COL ESEEESD SEL D 


GA 4O VD VD VD VD VD VI VM VD VD VM oo o o o CO 00 00 00 00 00 00 00 Y Y JI -2 12 12 -1 -2 -2 -2-20310101 010 O. O O O OUI UI Un Ul 


M" 


.6516e-01 
,15336-01 
. 8545e-01 
.9552e-01 
.0554e-01 
.1551e-01 
.2543e-01 
.3530e-01 
.4512e-01 
.5489e-01 
.6460e-01 
.7426e-01 
.8387e-01 
.9342e-01 
.0292e-01 
.1237e-01 
.2176e-01 
-31106501 
.4038e-01 
.4961e-01 
- 58/8801 
.6789e-01 
.7695e-01 
.8596e-01 
.9490e-01 
.0379e-01 
.1263e-01 
.2140e-01 
.3012e-01 
.3878e-01 
.4739e-01 
.5594e-01 
.6443e-01 
.7286e-01 
.8124e-01 
.8955e-01 
.9782e-01 
.0602e-01 
.1417e-01 
.2225e-01 
. 3029e-01 
. 3826e-01 
.4618e-01 
.5404e-01 
.6184e-01 
,6958e-01 
.7727e-01 
.8491e-01 
.9248e-01 
.0000e+00 


dd IUS SS] —J OO CO CO OO OO 00 00 00 00 OO OO CO CO CO CO CO CO CO CO CO CO OO OO CO CO CO CO OO CO «XO XO XO WOW WOO 


T 


T 


1 


.1985e-01 
.1694e-01 
.1399e-01 
.1101e-01 
.0799e-01 
.0493e-01 
.0184e-01 
.9872e-01 
.9557e-01 
.9238e-01 
.8916e-01 
.8592e-01 
. 8264e-01 
. /934e-01 
. 7601e-01 
.7266e-01 
.6928e-01 
.6587e-01 
.6245e-01 
.5900e-01 
.5552e-01 
.5203e-01 
.4852e-01 
.4499e-01 
.4144e-01 
.3787e-01 
.3429e-01 
.3069e-01 
.2708e-01 
.2305e6-DI 
.1981e-01 
.1616e-01 
.1249e-01 
.0882e-01 
.0513e-01 
.0143e-01 
.9773e-01 
.9401e-01 
.9029e-01 
.8657e-01 
.8283e-01 
. /909G-01 
. 7535e-01 
. 7160e-01 
.6785e-01 
.6410e-01 
.6034e-01 
.5658e-01 
.5282e-01 
.4906e-01 


a Gs LD UD UD UU Un LU LL LnLn Oh nn Y CYLL CY CY Cr y, ry O0 01 OY OY Oy -J -J - - - - - - ---- - --c0c0 


On 
ON 


fe 


P 


1 


.1201e-01 
.0563e-01 
.9919e-01 
.9270e-01 
.8616e-01 
.7958e-01 
.7297e-01 
.6631e-01 
.5963e-01 
.5291e-01 
.4616e-01 
.3939e-01 
.3260e-01 
25996 Om 
.1896e-01 
.1211e-01 
.0526e-01 
.9840e-01 
.9153e-01 
.8465e-01 
.7778e-01 
. 7090e-01 
. 6403e-01 
.5717e-01 
.5031e-01 
. 4346e-01 
. 3662e-01 
.2980e-01 
.2299e-01 
. 1620e-01 
.0944e-01 
.0269e-01 
.9596e-01 
.8927e-01 
EE RE 
75956- 01 
.6934e-01 
.6275e-01 
.5620e-01 
.4969e-01 
.4321e-01 
.3676e-01 
.3036e-01 
.2399e-01 
.1766e-01 
.1138e-01 
.0513e-01 
.9893e-01 
.9278e-01 
.8667e-01 


Ch Ch Ch Ch Ch Ch Ch Ch Ch Ch CN sl zl zl zl zl zl zl zl zl zl zl zl zl zl zl zl d zl d JO OO OO OO CO CO CO CO CO CO CO CO CO CO CO 


uh 
P; 


.8277e-01 
.7860e-01 
.7439e-01 
.7013e-01 
.6583e-01 
.6148e-01 
.5710e-01 
.5267e-01 
.4821e-01 
.4371e-01 
.3918e-01 
.3461e-01 
.3001e-01 
.2538e-01 
.2072e-01 
.1603e-01 
.1132e-01 
.9658e-01 
.0182e-01 
.9704e-01 
.9224e-01 
.8741e-01 
.8258e-01 
.7772e-01 
.7285e-01 
.6797e-01 
.6307e-01 
.5816e-01 
.5325e-01 
.4832e-01 
.4339e-01 
. 3845e-01 
.3350e-01 
.2855e-01 
.2360e-01 
.1865e-01 
.1370e-01 
.0874e-01 
.0379e-01 
.9884e-01 
.9390e-01 
.8896e-01 
.8402e-01 
. 7/909e-01 
.7417e-01 
.6926e-01 
.6435e-01 
.5946e-01 
. 5458e-01 
.4970e-01 


po po p p. pp. pp. p. pp. p. pp. p. B p B. p. p ES. p. p BH. HAHAHAHA jS. ES. BO B p. B). BS. BS. p. HEB. B. p. p. S. BH. [B p [D p p |D. |pb. |p pb pb nm] 


A 
Az 





.3023e+00 
.2853e+00 
.2692e+00 
.2538e+00 
.2392e+00 
.2253e+00 
.2120e+00 
.1994e+00 
.1873e+00 
.1759e+00 
.1649e-«00 
.1545e+00 
.1446e+00 
.1352e+00 
. 1262e+00 
.1176e+00 
.1095e+00 
.1018e+00 
.0944e+00 
.0874e+00 
.080 8e+00 
.0745e+00 
.0685e+00 
.0629e+00 
.0576e+00 
.0525e+00 
.0478e+00 
.0433e+00 
.0391e+00 
.0351e+00 
.0314e+00 
.0279e+00 
.0247e+00 
.0217e+00 
.0189e+00 
.0163e+00 
.0139e+00 
.0118e+00 
.0098e+00 
.0081e+00 
.0065e-«00 
.0051e+00 
.0039e+00 
.0028e+00 
.0019e-«00 
.0012e+00 
.0007e+00 
.0003e-«00 
.0001e+00 
.0000e-«00 


APPENDIX C. FANNO FLOW PROPERTIES 





















































EQUATIONS 
1. Dimensionless velocity 
n 1 
* p es 2 Y 2 EN 
dr Zip IE 5 vc) 
a Critical temperature ratio 
zw 
| llis yp 
n: y+ I 2 
3. Critical pressure ratio 
1 
TI 2 Ms ll >: 
M ml IT m|} 
4. Critical stagnation pressure ratio 
yrl 
Ee m E EE 2(Y-1) 
p: My+l 2 
5: Critical length ratio 
* 1-M? [y+1 Ii. 
ux] N EE In| M^ 1—— Je M | 
D yM 2Y Da | 2 














SE 


B. MATLAB CODE 


Farmo.m 
Fanno flow parameters 
Variables: 


M - Mach Number 
k - Gamma (Specific Heat Ratio) 
MSTAR - Dimensionless Velocity 
TRAT - Static Temperature Ratio 
PRAT - Static Pressure Ratio 
BIGPRAT - Stagnation Pressure Ratio 
FLDRAT - Critical Length Ratio 


Z ge oe ge o? oe oe oe dŷ d? o? o? o? o? o? 


[0:01 1 
k= 1.67; 
MSTAR - M.*(((2/(krl)) ^CEE((CCK-1)/2) ^M) ee 


TRAT = ( (27 (ice) * CF UO 2) MA SR IE N 


BIGPRAT = (1./M).*(((2/(k+1))*(1+(((k-1)/2).*M.^2))).^((k+1)/(2* ... 
AS 
FLDRAT = ((1-M.%2) ./(k.*(M.%°2) ) )+( (kt1)7 (2*k) ) - * Tog IER 


Y- [M; MSTAR ; TRAT ; PRAT; BIGPRAT; FLDRAT] ; 
fid = fopen('fanno.txt','w'); 


fprintf (fid, '%2.2f\t %4.4e\t %4.4e\t %4.4e\t %4.4e\n',Y); 
fclose(fid); 
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Table C.1. Fanno Flow (y=1.67) 


C) OO 0 OO CC9)€6)C c) O OO Q OO O OOOO OC OO OO oc o o 0 0 ' COOL OD00 


U1 U1 UI UI UI UI ll i ds UPS aS aS aS a Vë UJ UU CO C) UJ. UJ CJ Y WU PO BO BO PO CO PO PO FO CO ER. EP Pe i Un E. EE E. E iQ CO OY U1 i L2. BO. EB. (CO 





. 0000e+00 
.1554e-02 
.3107e-02 
.4657e-02 
.6204e- 02 
.7747e-02 
.9284e-02 
.0813e-02 
.2335e-02 
.0385e-01 
-15356 01 
.2684e-01 
.3832e-01 
.4978e-01 
.6123e-01 
.7266e-01 
.8408e-01 
.9548e-01 
.0686e-01 
.1821e-01 
.2955e-01 
.4087e-01 
.5216e-01 
.6342e-01 
.7466e-01 
.8588e-01 
.9706e-01 
.0822e-01 
. 1935e-01 
.3045e-01 
.4152e-01 
.5255e-01 
.6355e-01 
.7452e-01 
.8545e-01 
.9635e-01 
.0721e-01 
.1803e-01 
.2881e-01 
.3956e-01 
.5026e-01 
.6092e-01 
.7154e-01 
.8212e-01 
.9266e-01 
.0315e-01 
.1360e-01 
.2401e-01 
.3436e-01 
.4468e-01 
.5494e-01 


p) p p p p. pP p D P P p pP pP. p pP p pP pP AA AAA AAA AAA AA AAA pipmppp pmmppypmpmpbmnpmpmmmmpmbmgppmnd! 


T 
que 





.3350e-«00 
.3350e-«00 
.3348e+00 
.3346e400 
. 3343e-«00 
.3339e+00 
. 3334er00 
.3328e+00 
.3321e+00 
.3314e+00 
.3305e+00 
.3296e-«00 
.3286e+00 
.3275e+00 
.3263e+00 
.3250e+00 
.3236er00 
32226100 
.3207e+00 
.3190e+00 
.3173e+00 
.3156e+00 
.3137e+00 
.3118e+00 
.3097e+00 
.3076e+00 
.3054e+00 
.3032e+00 
.3008e+00 
.2984e+00 
.2959e+00 
.2934e+00 
.2907e+00 
.2880e+00 
.2852e+00 
.2824e+00 
.2795e-00 
.2765e+00 
.2734e+00 
.2703e+00 
.2671e-00 
.2638e+00 
.2605e+00 
.2571e-00 
.2537e+00 
.2502e+00 
.2466e+00 
.2430e+00 
.2393er00 
.2356e+00 
.2318e+00 


NO NO NO NO NO BO NO NO BO BO NO NO BO UJ UJ UJ UU UJ UJ UJ UJ CJ is GG 01 UI JO PR PPPRPPRN NUI 


dl 
NO 


P 
P* 





oo 


.1554e-«02 
.7767e-01 
.8508e-«01 
.8878e+01 
.3099e+01 
.9245e+01 
.6492e+01 
.4427e+01 
.2821e-«01 
.1535e+01 
.0483e+01 
.6054e+00 
. 8628e+00 
.2261e+00 
.6739e+00 
.1906e-«00 
.7639e-«00 
.3845e+00 
.0447e+00 
. 7388e+00 
.4618e-00 
.2099e+00 
.9796e+00 
. 7685e+00 
.5741e+00 
.3945e+00 
.2280e+00 
.0734e+00 
.9292e+00 
.7946e+00 
.6686e+00 
. 5503e+00 
.4391e+00 
.3344e+00 
.2355e+00 
.1420e+00 
.0535e+00 
.9696e+00 
.8899e+00 
.8141e+00 
.7420e-«00 
.6732e+00 
.6075e+00 
.5447e+00 
.4847e-00 
.4272e+00 
.3721e+00 
.3193e+00 
.2685e+00 
.2198e+00 


p p p p pP P pP. p. p. DD. pp. p IO. pO. jp. BÍ jj. pj. p. Bj ES: ND. BO BO. IO BO DO BO BO BO DO CU UJ CJ UO CJ iiS dS iS UT UI OY -J O0 x0 P OD. OD. fO Uf 


a 
P* 


[ 





OO 


.6235e+01 
.8123e+01 
.8755e+01 
.4073e+01 
.1265e+01 
.3943e+00 
.0593e+00 
.0589e+00 
.2817e+00 
.6607e+00 
.1533e+00 
.7310e+00 
.3744e+00 
.0692e+00 
.8052er00 
.5747er00 
.3718e+00 
.1919e+00 
.0313e+00 
.8871e-00 
.7571e-00 
.6392e+00 
.5319e+00 
.4339e+00 
.3440e+00 
.2614e+00 
.1852e+00 
.1147e-00 
.0494e+00 
.9886e-00 
.9321e+00 
.8794e-00 
.8301e+00 
. 7839er00 
.7406e+00 
. 7000e+00 
.6618e+00 
.6258e+00 
.5919e+00 
.5599e+00 
.5297e+00 
.5011e+00 
.4740e+00 
.4484e-00 
.4242e+00 
.4011e+00 
.3793er00 
.3585e+00 
. 3388er00 
.3201e+00 


DOMWWRPPRPRPPRPPPNONNNNWWWWRKRUUDA TWOP RPRPRPRPNNNWWHhRUNDWPRPNWOAEF UT 


AfL* 
D 





oo 


. 9803e+03 
. 4904e:03 
.5936e-«02 
.6874e+02 
. 3436er02 
.6147e+02 
. 1758e+02 
. 9155e+01 
.9706e-«01 
.5828e+01 
. 5588e+01 
. 7822e+01 
. 179 8e+01 
.7035e-«01 
.3206e+01 
. 0086e+01 
.7511e-«01 
.5363e+01 
.3555e+01 
.2018e+01 
.0704e+01 
.5705e+00 
.5879e+00 
. 7311e+00 
. 9800e+00 
.3186e+00 
.7335er00 
.2140e+00 
.7510e+00 
.3369e+00 
.9654e+00 
.6312e+00 
.3296e+00 
. 0569e+00 
.8097e+00 
.5852e-«00 
.3807e+00 
.1943e+00 
. 0239e+00 
. 8680e+00 
. 7251e+00 
.5939e+00 
.4733e+00 
.3623e+00 
.2601e+00 
.1658e+00 
.0787e+00 
.9823e-01 
.2379e-01 
.5488e-01 


Table C.1. Fanno Flow (y=1.67) (cont.) 


HD WWW ww WN VM VD VM VD OO OO OO OO OO OO 09 —J zl zl zl zl zl zs zl zl zl HAOOOOOOOOOOUY UIUYU UI 





.6516e-01 
. /533e-01 
.8545e-01 
.9552e-01 
.0554e-01 
.1551e-01 
.2543e-01 
. 3530e-01 
.4512e-01 
.5489e-01 
.6460e-01 
.7426e-01 
.8387e-01 
.9342e-01 
.0292e-01 
.1237e-01 
.2176e-01 
.3110e-01 
.4038e-01 
.4961e-01 
.5878e-01 
.6789e-01 
.7695e-01 
.8596e-01 
.9490e-01 
.0379e-01 
.1263e-01 
.2140e-01 
.3012e-01 
.3878e-01 
473986201 
.5594e-01 
.6443e-01 
.7286e-01 
.8124e-01 
.8955e-01 
.9782e-01 
.0602e-01 
.1417e-01 
.2225e-01 
.3029e-01 
. 3826e-01 
.4618e-01 
.5404e-01 
.6184e-01 
. 6958e-01 
, 7727e-01 
.8491e-01 
.9248e-01 
. 0000e+00 


PbebebperPRerEPREPREPPPPPPPPPPPPPPPEREPPPPPPRPPPPPPPPPPPPPPPPRPRPRPRP 


di 
Ts 





.2280e+00 
.2241e+00 
.2202e+00 
.2162e+00 
.2122e+00 
.2081e-«00 
.2040e+00 
.1998e+00 
.1956e+00 
.1913e+00 
.1870e+00 
.1827e+00 
.1783e+00 
.1739e+00 
.1695e+00 
. 1650e+00 
. 1605e+00 
. 1559e+00 
.1514e+00 
.1468e+00 
.1421e+00 
.1375e-00 
.1328e+00 
.1281e+00 
.1233e+00 
.1186e+00 
.1138e+00 
.1090e+00 
.1042e+00 
.0993e+00 
. 0944e+00 
. 0896e+00 
. 0847e+00 
.0798e+00 
.0748e+00 
.0699e+00 
.0650e+00 
.0600e+00 
05502100 
. 0501e-00 
.0451e+00 
.0401e+00 
.0351e+00 
.0301e+00 
.0251e-«00 
.0201e-«00 
.0151e+00 
.0100e+00 
. 0050e+00 
. 0000e+00 


PreobePrPREPRPPEPEPPRPPREPPPPPPPPPPEPPPREPRPPPPPRPRPPPPPPPPPPPREPRENNNDD 


ÓN 
cs 


P 
P" 





.1728e+00 
.1277e+00 
.0842e-«00 
. 0422e*r00 
.0018e+00 
.9627e+00 
.9250e-«00 
.8885e+00 
.8533e+00 
.8191e+00 
.7861e-00 
.7541er00 
. 7230er00 
.6929e+00 
.6637e+00 
.6354e+00 
.6078e+00 
.5811e+00 
.5551e00 
.5298e-00 
.5052e-«00 
.4813e-«00 
.4580e+00 
.4353e-«00 
.4132er00 
.3916e+00 
.3706e+00 
.3501e+00 
.3301e+00 
.3106e+00 
.2916e+00 
.2730e+00 
.2548e+00 
.2370er00 
.2197e+00 
.2028e400 
.1862e+00 
. 1700e+00 
. 1541e+00 
.1386e+00 
.1234e+00 
.1085e-«00 
.0940e-«00 
.0797e-00 
.0657e-«00 
.0521e+00 
.0387e+00 
. 0255er00 
.0126e+00 
.0000e+00 


PREP PP PPP PPP PPE PPP PPP RPRPRPRPRPPRPRPRPRPRPRPRPRPRPRPRPRPRPRPRPRPRPRPRPRPREFREF 


F 


I 


P; 





.3023e+00 
.2853e+00 
.2692e+00 
.2538e+00 
.2392e+00 
.2253e+00 
.2120er00 
.1994e+00 
.1873e+00 
.1759e+00 
.1649e+00 
.1545e-«00 
.1446e+00 
.1352e+00 
.1262e+00 
.1176e+00 
.1095e-«00 
.1018e-«00 
.0944e+00 
.0874e+00 
.0808e+00 
.0745e+00 
. 0685e+00 
.0629e+00 
.0576e+00 
.0525e+00 
.0478e+00 
.0433e+00 
.0391e+00 
.0351e+00 
.0314e+00 
.0279e+00 
.0247e+00 
.0217e+00 
.0189e+00 
.0163e+00 
.0139e+00 
.0118e+00 
.0098e+00 
.0081e+00 
.0065e«00 
.0051e+00 
.0039e+00 
. 0028e+00 
.0019e+00 
.0012e+00 
.0007e+00 
.0003e+00 
.0001e+00 
.0000e+00 


C XO UJ O0 H2 NJ) UJ U1 OY O0 IS EA kA kA RÄ AG GU Oz zs O0 X0 P2. ES i|). EB. [B|. ED. SS. BO. BO BO. BO UJ UJ UJ UJ 4» wm 4 UI Ul OY Ov - - 


AfL* 
D 





.9104e-01 
.3186e-01 
.7696e-01 
.2601e-01 
.7870e-01 
.3474e-01 
.9390e-01 
.5593e-01 
.2062e-01 
.8778e-01 
.5723e-01 
.2881e-01 
. 0237e-01 
.7778e-01 
.5490e-01 
.3361e-01 
.1383e-01 
,9543e-01 
.7834e-01 
.6246e-01 
.4772e-01 
.3405e- 01 
.2137e-01 
.0962e-01 
.8752e-02 
.8704e-02 
.9425e-02 
.0869e- 02 
.2992e-02 
.5752e-02 
.9111e-02 
. 3034e-02 
.7486e-02 
.2435e-02 
.7853e-02 
.3712e-02 
. 9986e-02 
.6651e-02 
. 3683e-02 
.1061e-02 
.7662e-03 
.7781e-03 
.0795e-03 
. 6535e-05 
.4843e-03 
.5572e-03 
.5802e-04 
.7362e-04 
.1530e-05 
. 0000e+00 


APPENDIX D. ISOTHERMAL FANNO FLOW PROPERTIES 























EQUATIONS 
1. Dimensionless velocity 
ep YM 
2 Critical stagnation pressure ratio 
ETES 
3 Critical stagnation temperature ratio 
laa iF d 
4 Critical length ratio 
we E + In (7?) 


B. MATLAB CODE 


iFanno.m 
Fanno flow parameters (isothermal flow) 
Variables: 


M - Mach Number 
k - Gamma (Specific Heat Ratio) 
PRAT - Static Pressure Ratio 
BIGTRAT - Stagnation Temperature Ratio 
BIGPRAT - Stagnation Pressure Ratio 
FLDRAT - Critical Length Ratio 


NO cP dO do? o? o9 o0 o9 o? o? o0 o0 o? 


B 


— LO. OL DUE 

MO Fr [./G 0 |. 

M - [Ml 1/sqrt(1.67) M2]; 

k= 1.67; 

PRAT = sqrt (k) .*™; 

BIGIRAT — (((2*5k)/ ((S*EK) O) CE (O EON Y DD 


BIGPRAT 


(1./PRAT).*((BIGIRAT).^(k/(k-1))); 

FLURAT - (1—(k.*(M. 2)))./(K.*(M. 2) leg PT ME 
Y- [M; PRAT; BTGTRAT ; BIGPRAT ; FLDRAT] ; 

fid - fopen('ifanno.txt','w'); 


forintf(fid, '%2.2f\t %4.4e\t %4.4e\t %4.4e\t %4.4e\n',Y); 
fclose(fid); 
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Table D.1. Isothermal Fanno Flow (y=1.67) 


ooo OO CG ee OO 2 OC O C91 CO 69 CO o (co CO COLCO 6 CC eee eo Oo eae E 


OY OY OY OY U1 Ul Ui Ul YU LI UI LI Va, VS Va S aS aS a as UJ UJ CJ CJ CJ CJ UJ NO NY INO IO DY DY PO BO AA kA kA kA kA kA kA Mä zl OU UJ ND P). o 





Ia PS 





View: Sie 


.0000e+00 
.2923e-02 
., 5846e- 02 
.8769e-02 
.1691e-02 
.4614e-02 
.7537e-02 
.0460e-02 
.0338e-01 
.1631e-01 
.2923e-01 
.4215e-01 
.5507e-01 
.6800e-01 
.8092e-01 
.9384e-01 
.0677e-01 
.1969e-01 
.3261e-01 
.4553e-01 
.5846e-01 
.7138e-01 
.8430e-01 
.9723e-01 
.1015e-01 
.2307e-01 
.3599e-01 
.4892e-01 
.6184e-01 
.7476e-01 
.8769e-01 
.0061e-01 
.1353e-01 
.2645e-01 
.3938e-01 
.5230e-01 
.6522e-01 
.7815e-01 
.9107e-01 
.0399e-01 
.1691e-01 
.2984e-01 
.4276e-01 
.5568e-01 
.6861e-01 
.8153e-01 
.9445e-01 
.0737e-01 
.2030e-01 
.3322e-01 
.4614e-01 


LO OO OO OO CO OO CO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OO OD OO OO OO OO OO OO OO OO OO OO OO OO CO CO CO CO CO CO CO CO CO CO 


P 


t 


P; 





.3292e-01 
.3295e-01 
.3303e-01 
.3317e-01 
.3336e-01 
.3362e-01 
.3392e-01 
.3428e-01 
.3470e-01 
.3518e-01 
.3571e-01 
.3629e-01 
.3694e-01 
-37036 01 
.3839e-01 
.3920e-01 
.4006e-01 
.4098e-01 
.4196e-01 
.4299e-01 
.4408e-01 
.4522e-01 
.4642e-01 
.4768e-01 
.4899e-01 
.5036e-01 
.5178e-01 
.5326e-01 
.5479e-01 
.5638e-01 
.5803e-01 
.5973e-01 
.6149e-01 
.6330e-01 
.6517e-01 
.6710e-01 
.6908e-01 
.7112e-01 
.7321e-01 
. /936e-01 
.7756e-01 
. 7/982e-01 
.8214e-01 
.8451e-01 
.8694e-01 
.8942e-01 
.9196e-01 
.9455e-01 
.9721e-01 
.9991e-01 
.0267e-01 
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L9 p p. DB p PS. p. p. B. BB. EA EA EA EA kA [O. [S). p. B. B. BS. (S. [EB E. E PO IO PO BO IO PO FO FO IOGUU C Yu y > ji U1 OY - p 00 iO E EB INO s 


T 


t 


T; 





co 


.9065e+01 
.4539e+01 
.6366e+01 
.2282e+01 
.8328e+00 
.2015e+00 
.0375e+00 
.1655e-«00 
.4882e+00 
.9472e+00 
.5053e+00 
.1378e+00 
.8274e+00 
.5620e+00 
.3326e+00 
.1323e+00 
.9561e+00 
.8000e+00 
.6607e+00 
.5358e+00 
.4232er00 
.3213e+00 
.2286er00 
.1440e+00 
.0665e+00 
.9953e+00 
.9297e+00 
.8692e+00 
.8131e+00 
.7611e-00 
.7127e-00 
.6676e+00 
.6256e+00 
.5863e+00 
.5496e+00 
.5151e+00 
.4828e+00 
.4524er00 
.4239er00 
. 3970er00 
.3717e-00 
.3478e+00 
.3253e+00 
.3041e+00 
.2840e+00 
.2651e+00 
.2472e-00 
.2302e+00 
.2142e+00 
.1990e+00 


LO UD Ch sl sl DV kA kA kA kA kA o EO BO ISONO NOLO C) CO 4S S U1 01 OY -J -2 OO OQ I B p p I IO FOCQU CO S U1 OY O0 lr EP ND OHUI 


AfL ' 
D 


oo 


.9783e+03 
.4887e-03 
.5784e+02 
.6733e+02 
.3304e+02 
60226402 
.1640e+02 
.8024e+01 
.8623e+01 
.4788e+01 
. 4586e+01 
. 6856e+01 
.0864e+01 
.6132e+01 
.2332e+01 
.9238e+01 
. 6689e+01 
.4565e+01 
.2179e+01 
.1264e+01 
.9698e+00 
.8565e+00 
. 8930e+00 
.0545e+00 
.3211e+00 
.6767e+00 
.1082e+00 
.6047e+00 
.1572e+00 
.7582e+00 
. 4015e+00 
. 0816e+00 
.7941e+00 
.5352e-«00 
.3014e+00 
. 0899e+00 
.8983e+00 
. 7245er00 
.5665e+00 
.4228e+00 
.2918e+00 
.1724e+00 
.0634e+00 
.6385e-01 
.7286e-01 
.8965e-01 
.1352e-01 
.4385e-01 
.8009e-01 
.2174e-01 


Table D.1. Isothermal Fanno Flow (y=1.67) (cont.) 
V PS 


mbocc coc ceoecccoe0 COCO Ce UO Oc co o oo CXCOC) CX GIC) CO C ESE ELE 


X 


pa pa p p) A kä kA kä kä kä kA kA kA kA kA kA kA kä kä kA kA kA kA ES: SQ «XO XO XO VM X0 40 30 00 00 00 00 00 00 00 O0 - - - -J -J -)-JOYO0Y0OY ON 





ur 


.5907e-01 
. /199e-01 
.8491e-01 
.9783e-01 
.1076e-01 
.2368e-01 
.3660e-01 
.4953e-01 
.6245e-01 
.7537e-01 
.8829e-01 
.0122e-01 
.1414e-01 
.2706e-01 
.3999e-01 
.5291e-01 
.6583e-01 
.7875e-01 
.9168e-01 
.0460e-01 
.1752e-01 
. 3045e-01 
.4337e-01 
.5629e-01 
.6921e-01 
.8214e-01 
.9506e-01 
. 0000e+00 
.0080e+00 
.0209e+00 
.0338e+00 
.0468e+00 
.0597e+00 
.0726e+00 
. 0855e+00 
. 0984e+00 
.1114e+00 
.1243e+00 
.1372e+00 
.1501e+00 
.1631e+00 
.1760e+00 
.1889e+00 
.2018e+00 
.2147e+00 
.2277e-00 
.2406e+00 
.2535e+00 
.2664e+00 
.2794e+00 
.2923e*r00 


p p p) p p pm|B HB mp pmP pP P P P P P P H PH HP H ddr Mä VD VD XO XO XO XO XO XO XO UO XO AO XO XO LO XO KO LO LO LO KO LO kO LO LO 


P, 
P; 





.0549e-01 
.0837e-01 
.1130e-01 
.1428e-01 
.1732e-01 
.2042e-01 
.2357e-01 
.2678e-01 
.3005e-01 
.3337e-01 
.3674e-01 
.4018e-01 
. 4366e-01 
.4721e-01 
. 5081e-01 
.5446e-01 
.5817e-01 
.6194e-01 
.6576e-01 
.6964e-01 
.7358e-01 
.7757e-01 
.8161e-01 
.8571e-01 
. 8987e-01 
. 9408e-01 
. 9835e-01 
. 0000e+00 
.0027e+00 
.0071e+00 
.0115e+00 
. 0160e+00 
. 0205e+00 
. 0251e+00 
. 0298e+00 
.0345e+00 
.0393e+00 
.0441e-«00 
. 0490e+00 
. 0539e+00 
. 0589e+00 
. 0640e+00 
.0691e+00 
.0742e+00 
.0795e+00 
.0847e+00 
.0901e+00 
.0955e+00 
.1009e+00 
.1064e-«00 
.1119e+00 
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LA Lä ka ka XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO VD XO XO ID. I2. [S. [) [O. EO. E. i. BB. p. p. EB. p. p. ES ES. p. BÍB. BÍ. EB. [| E]. EO. E. i. E. EB. ES 


T 
T; 





.1847e+00 
.1711e+00 
.1583e+00 
.1461e+00 
.1347e-00 
.1238e+00 
.1135e-«00 
.1038e+00 
.0947e+00 
.0860e+00 
.0779e+00 
.0702e+00 
.0630e+00 
.0562e+00 
. 0498e+00 
.0439e+00 
.0383er00 
.0331e+00 
.0282e+00 
.0237e+00 
.0195e+00 
.0157e+00 
.0121e+00 
.0089e+00 
.0059e+00 
.0032e+00 
. 0008e+00 
. 0000e+00 
.9872e-01 
.9685e-01 
.9523e-01 
.9386e-01 
.9273e-01 
.9184e-01 
.9117e-01 
.9073e-01 
.9050e-01 
.9049e-01 
.9069e-01 
.9110e-01 
.9170e-01 
.9251e-01 
.9350€e DI 
.9469e-01 
.9607e-01 
.9763e-01 
.9937e-01 
.0013e+00 
.0034e+00 
.0057e+00 
.0081e+00 


BPrRIDMOTWAAHAUABWWNHONKPFPPUWODPBNDWRFPORAMEPARAIRPRPNNWHEUNADAUWORPRPRPERNNNNWWEHEE 


AfL* 
D 





.6834e-01 
.1948e-01 
.7480e-01 
.3396e-01 
.9666e-01 
.6263e-01 
.3162e-01 
.0340e-01 
.7776e-01 
.5451e-01 
.3348e-01 
.1451e-01 
.7451e-02 
.2169e-02 
.8542e-02 
.6454e-02 
. 5801e- 02 
.6487e-02 
. 8420e- 02 
.15196202 
.5707e-02 
.0913e-02 
.0696e-03 
.1167e-03 
.9971e-03 
.5768e-04 
.9226e-05 
.0000e+00 
.2575e- 04 
.4442e-04 
.1653e-03 
.0510e-03 
.4667e-03 
.3797e-03 
.2760e-02 
.6578e-02 
.0808e-02 
.5424e-02 
.0404e- 02 
.5724e-02 
. 1365e-02 
.7306e-02 
. 3530e-02 
.0019e-02 
.6758e-02 
.3730e-02 
.0922e-02 
.8320e-02 
. D911& 305 
. 0368e-01 
.1163e-01 


APPENDIX E. RAYLEIGH FLOW PROPERTIES 























EQUATIONS 
1. Dimensionless velocity 
* (yr1)M^ 
we B Un 
p 1+ yM 
2: Critical temperature ratio 
T (y+1)M 
T' (1+ y 
3. Critical pressure ratio 
P yi 
P* 1a yM? 
4. Critical stagnation pressure ratio 
zm El 2 EG 
P* 1ryM (YT! 2 
F Critical stagnation temperature ratio 


T, 2(y+1\M’ ET 
IU earl 
Tr (xa 
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B. MATLAB CODE 


Rayleigh.m 
Rayleigh flow parameters 
Variables: 


M - Mach Number 
k - Gamma (Specific Heat Ratio) 
MSTAR - Dimensionless Velocity 
TRAT - Static Temperature Ratio 
PRAT - Static Pressure Ratio 
BIGPRAT - Stagnation Pressure Ratio 
BIGIRAT - Stagnation Temperature Ratio 


r= oP a oP oP oP dO oO dO oP oe? o9 o9 oe oe 


[0:01 
k= 167: 
MSTAR = ((k+1).*(M.%2)) ./(1+(k.*(M.%2))); 
TRAT- (en) 2). "M. 2). de ds Ma SE ET 
PRAT = Kl) ML Mo UT SUME 
BIGPRAT = (PRAT) -*(((2/ (k+l) )* Ce E IU 
BIGTRAT = 


((2*(kaed).*(M.2)) ./ (OER 2 AA AE eee 
(os dIE M 2) 


Y= [M; MSTAR ; TRAT; PRAT; BIGPRAT; BIGIRAT] ; 
fid = fopen( 'rayleigh.txt', 'w'); 


fprintf (fid, '$2.2fNt %4.4e\t %4.4e\t %4.4e\t $4.4eNt %4.4e\n',Y); 
fclose(fid); 
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Table E.1. Rayleigh Flow 


CC OO0O0OOcUOCuoDooo0ooooodcocdocodguodocododdodgoododccododoodocdocoododoccocdotddooddoccd 


A 


DS ës ds a RB Gs Y UU. CU. UJ. UJ. UJ GJ NO BO NO NO BO BO BO BO ES. E. ES. F2. F2 ES ES HRS ES ES XO COO -J OY Ul Ul d» UJ C) BPO PO ES F3 «X0 Oh P2 ES. P2. COD 





. 0000e+00 
.6696e-04 
.0673e-03 
.3994e-03 
.2606e-03 
64 /2e- 05 
. 5546e-03 
.2977e-02 
.6907e-02 
.1338e-02 
.6261e-02 
.1667e-02 
.7545e-02 
. 3884e- 02 
.0673e-02 
.7899e-02 
.5550e-02 
. 3610e-02 
.2067e-02 
. 0907e-02 
.0011e-01 
.0967e-01 
.1956e-01 
.2978e-01 
.4030e-01 
.5110e-01 
.6218e-01 
273529] 
.8509e-01 
.9689e-01 
. 0890e-01 
.2110e-01 
.3348e-01 
.4602e-01 
.5871e-01 
.7153e-01 
.8446e-01 
.9751e-01 
.1064e-01 
.2385e-01 
.3712e-01 
.5045e-01 
.6381e-01 
.7721e-01 
.9062e-01 
.0404e-01 
.1746e-01 
.3086e-01 
.4424e-01 
.5759e-01 
.7090e-01 


CO OO OO CO CO ERR el sl «J «J «J OY OY Ov Oy OY Ui Ul Ul Ui i d» a UD Y uy ud NN NDNDNNDN EA EB. E E, E i: CO QY U1 i LO FO E E. OY P2. O 





. 0000e+00 
.1265e-04 
.8478e-03 
.3968e-03 
.1346e-02 
.7674e-02 
.5358e-02 
.4367e-02 
.4665e-02 
.6213e-02 
.8966e-02 
.2876e-02 
. 7/891e-02 
.1396e-01 
.3101e-01 
.4899e-01 
.6784e-01 
.8749e-01 
.0787e-01 
.2892e-01 
.5056e-01 
.7273e-01 
.9536e-01 
.1838e-01 
.4172e-01 
.6532e-01 
.8910e-01 
.1301e-01 
.3699e-01 
.6097e-01 
.8489e-01 
.0870e-01 
.3236e-01 
. 5580e-01 
.7898e-01 
.0185e-01 
.2438e-01 
.4653e-01 
. 6826e-01 
159536-DT 
.1032e-01 
.3060e-01 
.5034e-01 
.6953e-01 
.8814e-01 
.0616e-01 
.2358e-01 
.4037e-01 
,5655e-01 
, /209e-01 
.8699e-01 


GA L2 p2 ES LAS p2 BO BO BO BO BO BO BO NO BJ BO BO BO BO BO BO BO DJ bO bJO BJ DO BJ BJ DO bO BO bO BO BO BO BO BO BO BO BO BO BPO BO PO PO BO BO BO BO BO 


ON 
N 


P 
P* 





.6700e+00 
. 6696e*r00 
.6682e+00 
.6660e+00 
.6629e«00 
.6589e+00 
. 6540e+00 
.6483e+00 
.6418e+00 
.6344e+00 
.6261e+00 
.6171e+00 
.6073e+00 
.5967e+00 
. 5854e+00 
.5733e+00 
. 5605e+00 
.5471e+00 
.5329e+00 
.5182e+00 
.5028e+00 
.4869e-«00 
.4703e«00 
.4533e+00 
.4357e+00 
.4177e+00 
.3992e+00 
.3802e+00 
.3609e+00 
.3412e+00 
.3211e+00 
.3008e+00 
.2801e+00 
.2591e+00 
.2380e+00 
.2165e+00 
.1949e+00 
.1732e+00 
.1512e+00 
.1292e+00 
.1070e+00 
. 0848e+00 
.0624e+00 
.0401e+00 
.0177e+00 
.9953e+00 
.9729e*00 
.9505e-«00 
.9281e+00 
.9058e+00 
.8836e+00 


PPP PE PEP S. B. E. ES. EB. ES. F2. ES. F2. S. ES. ES. PRP RP RP PRP PRP F2. 2. ES. RP PRP ES. RP RP RP RP RP RP ES. RP RP RP RP E2. E2. F2. E. E2 


R 


t 


P; 





.2994e+00 
.2993e+00 
.2990e+00 
.2984e+00 
.2977e-00 
.2967e+00 
.2955er00 
.2941e+00 
.2925e+00 
.2908e+00 
.2888e+00 
.2866e+00 
.2842e+00 
.2816e+00 
.2789e-00 
.2760er00 
.2729e-00 
.2697e+00 
.2663e+00 
.2628e+00 
.2591e+00 
.2553e+00 
.2514e+00 
.2474e-00 
.2432e+00 
.2390e+00 
.2346er00 
.2302e+00 
.2257e+00 
.2211e+00 
.2164e+00 
.2117e+00 
.2070e+00 
.2022e+00 
.1973e+00 
.1925e+00 
.1876e+00 
.1827e+00 
.1778e+00 
.1728e+00 
.1679e+00 
.1630e+00 
.1581e+00 
.1533e+00 
.1484e+00 
.1436e+00 
.1388e+00 
.1341e+00 
.1294e+00 
.1248e+00 
.1202e+00 


-] «JOY OY OY OY OY OY Un Ul G1 Ul Un Un i iS dS ds se Y WWW WU YWYNNNNNNDNR?R E E E EB EB :Q OO -2ONU1 42 C0 lO F2 E CO S lO U1 (D 


T 


t 


T; 





. 0000e+00 
. 3384e-04 
.1334e-03 
. 7/930e-03 
.5031e-03 
:3250e8202 
.9018e-02 
.5785e-02 
. 3529e- 02 
.2221e-02 
. 1833e-02 
.2331e-02 
.3681e-02 
.5843e-02 
.8779e-02 
.1245e-01 
.2680e-01 
.4180e-01 
.5740e-01 
.7355e-01 
.9020e-01 
.0731e-01 
.2483e-01 
.4271e-01 
.6091e-01 
. 7937e-01 
.9806e-01 
.1693e-01 
.3593e-01 
.5502e-01 
.7416e-01 
.9332e-01 
. 1245e-01 
. 3152e-01 
.5049e-01 
69336201 
.8801e-01 
.0650e-01 
.2478e-01 
.4282e-01 
.6059e-01 
.7808e-01 
. 9527e-01 
.1213e-01 
.2866e-01 
.4483e-01 
.6064e-01 
.7608e-01 
.9113e-01 
. 0579e-01 
.2005e-01 


Table E.1. Rayleigh Flow 


x 


GA VD VD AO AO AO AO AO XO XQ XQ iO iO O0 000000000000000000 12-12 -)2 -2 -2 -2 -2 -2 -JOYOYOY OY OY OYOY OY OY U1 U1 C101 Ur Ur Un i e 


M'- 


p 
p 





.8416e-01 
.9737e-01 
. 1052e-01 
.2360e-01 
.3660e-01 
.4952e-01 
.6236e-01 
.7510e-01 
.8775e-01 
. 0030e-01 
.1274e-01 
.2508e-01 
.3730e-01 
.4941e-01 
.6140e-01 
.7328e-01 
.8503e-01 
.9665e-01 
. 0O815e-01 
.1952e-01 
.3076e-01 
.4187e-01 
.5285e-01 
.6370e-01 
.7441e-01 
.8499e-01 
.9544e-01 
.0576e-01 
.1594e-01 
.2599e-01 
. 3590e-01 
.4568e-01 
.5533e-01 
.6485e-01 
.7424e-01 
8350 OU 
.9262e-01 
.0162e-01 
.1050e-01 
.1924e-01 
.2786e-01 
.3636e-01 
.4473e-01 
: 92986-01 
.6111e-01 
.6912e-01 
.77026-01 
.8479e-01 
.9245&e-01 
.0000e+00 


pa pa p2 pp pa. pa. ES. pp. pp. pp. pp. po. p S. pO. p. B. RS. S. ES. ES. S. EO. ||. B. ES. B. B. ES. S. DO. ES. ES. ES. ES. ES. S. ES. EB. EB. ED: NO XO XO XO XO XO O AO AD 


(y=1.67) (cont.) 


n 
que 





.0124e-01 
.1486e-01 
.2783e-01 
.4017e-01 
.5186e-01 
. 6293e-01 
. /336e-01 
.8318e-01 
.9239e-01 
.0010e+00 
.0090e+00 
.0165e+00 
.0233e+00 
.0296e+00 
.0354e+00 
.0406e+00 
.0454e+00 
.0496e+00 
.0533e+00 
. 0565e+00 
.0593e+00 
.0617e+00 
. 0636e+00 
.0651e+00 
.0662e+00 
.0669e+00 
.0672e+00 
.0671e+00 
.0667e+00 
.0660e+00 
.0650e+00 
.0636e+00 
.0620e+00 
.0600e+00 
.0578e+00 
.0554e+00 
.0527e+00 
. 0497e+00 
.0466e+00 
.0432e+00 
.0396e+00 
.0359e+00 
.0319e+00 
.0278e+00 
.0235e+00 
.0191e+00 
.0145e+00 
. 0098e+00 
. 0050e+00 
. 0000e+00 


PrEbEPPPRPPPEPRPRPPPPREPPPPPPEREPRREPPPRPPREPPPPPPPPRPEPREPPREPEPPPPRPRER 


ON 
de 


le 
P* 





.8614e+00 
.8394e+00 
.8174e+00 
.7956e«00 
.7739e-00 
.7523e400 
.7309e+00 
. 7096e+00 
.6885e+00 
.6675e+00 
.6467e+00 
.6261e+00 
.6057e+00 
.5855e+00 
.5655e+00 
. 5456e+00 
.5260e+00 
.5066e+00 
.4874e+00 
.4684e+00 
.4496e-00 
.4311e+00 
.4127e+00 
.3946e+00 
.3767e+00 
.3591e+00 
.3416e+00 
.3244e+00 
.3074e+00 
.2906e+00 
.2740e+00 
.2577e-00 
.2416e+00 
.2257e+00 
.2100e+00 
.1946e+00 
.1793e+00 
.1643e+00 
.1495e+00 
.1349e+00 
.1205e+00 
.1063e+00 
.0923e+00 
.0785e+00 
.0649e+00 
.0516e+00 
.0384e+00 
.0254e+00 
.0126e+00 
. 0000e+00 


PEPbeobobRbPerPRPPRPPPPPPRPPPPPREPPPPERPRPRPPPPPPPPPPPPPPPPPRPRRRRPR 


d 


t 


P; 





.1156e+00 
.1112e+00 
.1067e-00 
.1024e+00 
.0981e-«00 
.0939e+00 
.0898e+00 
.0857e+00 
.0817e+00 
.0778e+00 
.0740e+00 
.0703e+00 
.0667e+00 
.0631e+00 
.0597e+00 
.0563e+00 
.0530e+00 
.0499e+00 
.0468e+00 
.0438e+00 
.0409e+00 
.0381e+00 
.0354e+00 
.0328e+00 
.0303e+00 
.0279e+00 
.0256e+00 
.0234e+00 
.0213e+00 
.0193e+00 
.0174e+00 
.0156e+00 
.0139e+00 
.0123e+00 
.0108e+00 
. 0094e+00 
.0081e+00 
.0069e+00 
.0058e+00 
.0048e+00 
. 00396 UO 
.0030e+00 
.0023e+00 
.0017e+00 
.0012e+00 
.0008e+00 
.0004e+00 
.0002e+00 
.0000e+00 
. 0000e+00 


GA VD XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO XO 00 00 00 Q0 O0 O0 00 00 00 00 JJ Y A -J 


T 


í 


T; 





.3391e-01 
.4736e-01 
.6041e-01 
.7304e-01 
.8526e-01 
.9707e-01 
.0847e-01 
.1946e-01 
.3005e-01 
.4024e-01 
.5003e-01 
.5943e-01 
.6845e-01 
.7709e-01 
.8535e-01 
.9325e-01 
. 0079e-01 
.0798e-01 
.1482e-01 
.2133e-01 
.2751e-01 
.3337e-01 
.3892e-01 
.4416e-01 
.4911e-01 
.5377e-01 
.5816e-01 
.6227e-01 
.6612e-01 
.6972e-01 
.7307e-01 
.7619e-01 
.7907e-01 
.8174e-01 
.8418e-01 
.8643e-01 
.8847e-01 
.9032e-01 
.9199e-01 
.9348e-01 
.9480e-01 
.9595e-01 
.9695e-01 
.9779e-01 
.9849e-01 
.9905e-01 
.9947e-01 
. 9977e-01 
.9994e-01 
.0000e+00 


APPENDIX F. FLOW CALCULATION CODE 


flowcalc.m 


A program to estimate the theoretical performance of a thermothrottle 
as a function of temperature for various values of throttle length 
and diameter and inlet pressure. 


oO o0 o9 o9 o? o0 oC 


Initial Conditions 
WO = 275. 

Mi. 1EO: 

eel rt) - 2/5:10:37b; 
PtO - 10000; 


Y 

| 
i 
On 
“dg 


ac 
2003; 

Pal (Pa) 4, 
MEE: 

Gon 525; 


1 
0 


Isentropic Flow to Station 1 


cere DN PORBRE 


E 


— Pt0/(R*TCO); 

Vstar —- sgrt(((2*k)/(krl))*R*TtO); 
rho0star = rho0* ( (2/(k+1) )*(1/(k-1))); 
mdotstari = rho0star*A*Vstar; 


Fanno Flow to Station 2 


Fh œ X oP 


IAS pl) 355) *(2-146/mgDESEar1) ; 

yT = (4*£* (L1+L2+L3) ) /D; 

ML = sart((1+(yT/2)- (sqrt ( (25*yT^2) - (100*yT) 410) /10) ) /k) ; 

pde cho (EE (dd) AE ML AD (k/ (k-1) )); 

rholstar - rhol*scrt(k)*Ml; 

mdotstarf - rholstar*A*Vstar; 

BE E: 

sE = yT ((4*E1*L1)/D); 

M2 sqrt ( (1- (y1/2)- (sqrt ( (25*y1^2) - (100*y1) 4-10) /10) ) /k) ; 

PEI at (kk M2) ) ~(((((2*k) /((3*k)-1))*... 
)) 
(KS 


REES ET ERC EC Mp... 
OSOD (Ta (M1) /2) XML 2) )) ^(Kk/(K—L))))5 


Rayleigh Flow to Station 3 


of dP ye 


69 


TtRatiol — (TES FO N 
((2* (k+1) * (NET AI 
TtRatiol; 
sqrt( ( (-sart (- ((k^2) *z1) - (2*k*z1)- (k^2) - (2*k) -1-z1) ) ... 
—(k*z1)-k-1)/(((xX^2)*z1)-(k^2)41)); 
Pt3 = Ee Ke EE 
(17 (((k-1) /2) *M3 2) ee APTA AAA 
(C7 Occ) y^ CL (UK MA i TUE 
rig) a DESCH -; 
rho3star = rho3*( ( (k+1) * (MB%2) ) / (1+ (k* (MB%2) ))); 
mdot3star - rho3star*A*Vstar; 


zi 
M3 


Fanno Flow to Station 4 


( (48*p1) /1E8) * ( (2 .146+ ( .00714* (Tt3-275) ) ) Amdot3star) ; 
(d*E3 ELE 


Station 4 to Station 1 (and back) 


dë dë oe KK rh ce oe o? 
LA LA 
WM 


ba. ES 

y3a = y3; 

Ma = sart( (1+ (y3a/2)- (sgrt ((25*y3a^2)+(100*y3a)+10)/10))/k); 

TtRatio2 = (Tt2/Tt3) * ((2* (k+1) * (MBa%2) ) / ( (1+ (k* (MBa%2) ))%2)*... 
(1+ ( ( (k-1) /2) *MBa%2))); 

z2 = TtRatio2; 


M2a = sart(((-sqrt (-( (k*2) *z2) - (2*k*z2)+ (k*2)+(2*k)+1—-z2))... 
-(k*z2 Fkri) / ( (0 2) 22) (KP) EDE 

fla = ((48*pi)/1E8)*(2.146/mdot3star); 

yla = log(k*(M2a^2))—((l-(k*(M2a^2)))/(k*(M2a^2)))*((4*fla*Ll)/D); 

Mia - scrt((l-(yla/2)-(sgrt((25*yla^2) K(l00 la DTE 


Pt2a = Ptl*((d/ (sort (k) Ma) TUE = oe 
Lk OKE1) /2) Ma2) (ef (ER (Scudo NN 
(CCCU2KE/ CC3*K) 1) ) *CEECCCOK-D)72]8 Mae 2E LIE 

Pt3a = Pt2a*(((k+1) /(1+(k* (MBa%2) ))) *(((2/ (k4+1))*... 
(TA) *M3a^2) ) ) ^ (K/ (k-1) AU de dik a 2 ee 
(((2/ (Kk+1) )* (1+ (((k-1) /2) *M2a%2) ) )*(k/ (k-1)))); 

rho3a = Pt3a/(R*Tt3); 

rhodastar = rho3a* ( ( (k+1) * (MBa%2) ) / (1+ (k* (MBa%2)))); 

% 

% Theoretical mass flow rate 

% 


mdot3astar = rho3astar*A*Vstar 


flb = ((48*p1)/1E8) * (2 .146/moet3 star) - 

ylb - (4*flb*Ll)/D; 

f3b = ((48*p1i) /1E8) * ((2.146+( .00714* (TE3—275) ) ) /mdot 3astar |; 
YD = (AE ES NE: 

end 
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